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REFINED SCATTERING AND HERMITIAN SPECTRAL THEORY 
FOR LINEAR HIGHER-ORDER SCHRODINGER EQUATIONS 

V.A. GALAKTIONOV AND I.V. KAMOTSKI 

Abstract. The Cauchy problem for a linear 2TOth-order Schrodinger equation 

P^ I (0.1) Mt = — i (-A)''"u in M^ x R+, u|t=o = "o; m>l is an integer, 

is studied, for initial data mq in the weighted space L^. (M^), with p*{x) — e'^l" and 
^ ' a = 2^73]- G (1,2]. The following five problems are studied: 

(I) A sharp asymptotic behaviour of solutions as t ^ +oo is governed by a discrete 
spectrum and a countable set $ of the eigenfunctions of the linear rescaled operator 

B = -i (-A)" + J-y.V + ^J, with the spectrum a(B) = {A^ = -M, |/3| > O}. 

Ph \ (II) Finite-time blow-up local structures of nodal sets of solutions as i — >■ 0~ and forma- 

■^^ ■ tion of "multiple zeros" are described by the eigenfunctions being generalized Hermite 

pH I polynomials, of the "adjoint" operator 

^ ; B* = -i (-A)" - ^y-y, with the same spectrum cr(B*) = cr(B). 

Applications of these spectral results also include: (III) a unique continuation theo- 
rem, and (IV) boundary characteristic point regularity issues. 

Some applications are discussed for more general linear PDEs and for the nonlinear 
^ ' Schrodinger equations in the focusing ("-I-") and defocusing ("— ") cases 

f^ ■ ut = -ii-A)"'u±i\u\P-^u in M^ x M+, where p > 1, 

^D ' as well as for (V) the quasilinear Schrodinger equation of a "porous medium type" 

ut = -i{-A)"'{\u\''u) in R^ x R+, where n > 0. 

f^ ' For the latter one, the main idea towards countable families of nonlinear eigenfunctions 

is to perform a homotopic path n -^ 0+ and to use spectral theory of the pair {B, B*}. 
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1. Introduction: duality of global and blow-up scalings, Hermitian 
spectral theory, and refined scattering 

1.1. Basic Shrodinger equations and key references. Consider the linear 2mth- 
order Schrodinger equation (the LSE-2m), with any integer m > 1, 

(1.1) ut = -i(-A)'"n in M^ x M+, u\t=o = mq, 
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where A is the Laplace operator in M^, for initial data Uq in some weighted L^-space, to 
be introduced. Here m = 1 corresponds to the classic Schrodinger equation 

(1.2) iut = -Au in M^ x M+, 

which very actively entered general PDE theory from Quantum Mechanics in the 1920s. 

It is not an exaggeration to say that, nowadays, linear and nonlinear Schrodinger type 
equations are the most popular PDE models of modern mathematics among other types 
of equations. In Appendix A, we present a "mathematical evidence" for that by using 
simple data from the MathSciNet. It not possible to express how deep is mathematical 
theory developed for models such as fll.21) . fll.ip . and related semihnear ones. We refer to 
well-known monographs [551 ITU] , which cover classes of both linear and nonlinear PDEs. 

Concerning the results that are more closely related to the subject of this paper, we 
note that scattering L^ and L^'^^ theory for (11.21) has been fully developed in the works by 
Stein, Tomas, Segal, Strichartz in the 1970s with later further involved estimates in more 
general spaces by Ginibre and Velo, Yajima, Cazenave and Weissler, Montgomery-Smith, 
Keel, Tao, and many others; see ^J and [91j for references concerning these, as well as 
optimal L^'^ estimates for the non-homogeneous equation 

(1.3) iut = -Au + F{x,t) in M^ x M, 

as well as more recent papers [5Ql |82], |99] . 

The 27Tith-order counterpart (II. ip was also under scrutiny for a long period. We refer to 
Ablowitz-Segur's monograph [1], Ivano-Kosevich [13], Turitsyn [89], Karpman [H], and 
Karpman-Shagalov [35] for physical, symmetry, and other backgrounds of higher-order 
Schrodinger-type semilinear models (see also [100] for extra motivations from nonlinear 
optics), [74] for first existence and uniqueness results, and more recent papers [3 El [T31 
E31 ESI ESI ESI ESI ESI IIOI] as an account for further applied and rigorous research, as 
well as other earlier key references and survey, in this fundamental area of modern PDE 
theory. 

1.2. Discrete real spectra, "Hermitian spectral history", and our motivation. 

Actually, the developed here refined scattering theory is rather general, so our methods 
can be applied not only to the Schrodinger equations such as (II. ip . but also to practically 
any linear evolution PDEs with constant or perturbed smooth coefficients and classic 
solutions. Formulating the approach rather loosely, we claim that, for (II. ip . the most 
principal part is played by the spectral theory for the following rescaled operator: 

(1.4) B* = -i{-Ayr-^yVy, where y = j^f^ for t<0 

is the corresponding Sturmian blow-up backward spatial variable at the focusing point 
(x,t) = (0,0~). In 1836, C. Sturm used the backward variable y = -^, t — )• 0~, for the 
heat equation with a potential, 

(1.5) ut = u^x + q{x)u, 

in his seminal paper [85] [The pioneering work [SS] was practically fully forgotten for 
about 150 yearsl-practically until the 1980s,-plausibly the most amazing and striking 
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such an example in the whole history of mathematics ever] , where he formulated his two 
fundamental theorems on zeros sets of solutions u = u{x,t) of fll.Sp . This remarkable 
history, with many extensions about, can be found in [23l, Ch. 1] with precise statements 
of Sturm's results of the 1836 written in his original notations. 

The operator (11.41) and its adjoint below are then respectively defined in weighted L^- 
spaces, with a special "radiation conditions at infinity" to be specified. More than half of 
the paper is devoted to the study of (11.41) and its "adjoint" operator 

(1.6) B = -ii-Ayr + ^yVy + £l, y=^, t>0 {the forward vanable). 

Indeed, (II. 4p is a perturbation of the original one in (II. ip . Though the perturbation is 
of the first order, the coefficient y therein is unbounded as y — )■ +oo, so this changes the 
natural space L^(M^) and moves the operator into an essentially weighted metric. An 
amazing property of (11.40 is that, being properly defined, it has the discrete spectrum 

(1.7) (t(B*) = {A/3 = -||, |/3| =0,1,2, ...} (/3isamultiindexinM^), 

and all the eigenfunctions {'0«(y)} are finite polynomials {generalized Hermite ones). 
These properties directly match the classic results for the heat equation for m = 1: 

ut = u,, =^ B* = Dl - lyDy, with (x{B*) = {Xp = ~^, I >0}, 
so(B*)* = B* in L%{R), p*{y)=e-4, B = Dl + lyDy+'^I, 

where B is defined in the adjoint space L% with p = 4r, and (B)* = B* in the dual 

metric of L^, etc. In modern language, for (II. Sp . the spectrum (II. 7p with m = 1 and 
the classic Hermite polynomials (introduced in detail about 1870) as eigenfunctions were 
already constructed by C. Sturm in 1836 [85]. As we mentioned already, this led Strum 
to formulate his two fundamental theorems on the structure of multiple zeros of solutions 
of parabolic equations and on nonincrease in time of the zero number (or sign changes 
of solutions); we refer again to [23], Ch. 1] for a full history and key further references 
and extensions. The operator B* in (II. 8p . admitting a natural A^- dimensional extension 
similar to (II. 4p . remains one of the key objectives in general theory of linear self-adjoint 
operators; see Birman-Solomjak's monograph [9j. 

The spectral results for (11.80 and their consequences for the asymptotic behaviour for 
second-order parabolic equations are classic and well-known since the 1830s, with further 
extensions as orthonormal polynomial families by Hermite himself from the 1870s. 

However, and this looks like a truly amazing fact, a direct extension of such classic 
results to other classes of PDEs took a lot of time. For instance, similar spectral theory 
for the ID hi-harmonic equation (cf. with the one in the first line in (II. 8p ) 

(1.9) ut = -u^^^^ =^ B* = -Dl-\ yDy, with a{B*) = [X^ = -|, / = 0, 1, 2, ...}, 

etc., was developed in 2004 [1^, i.e., 168 years later after Sturm's pioneering discovery 
for m = 1 in 1836! As a certain (but seems not that convincing) excuse, note that the 
operator B* in (II. 9p is not self-adjoint in no weighted space, though keeps having discrete 
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real spectrum, polynomial eigenfunctions (naturally called generalized Hermite ones), and 
a number of other nice and typical from self-adjoint theory properties. 

Therefore, our goal is to show that similar issues remain true for our rescaled Schrodinger 
operators f ll.4p and f ll.6p (so we are talking about a spectral pair {B,B*} of non-self- 
adjoint ones), which clearly have analogous structures, though the mathematics becomes 
essentially more involved than for (II. 9p . to say nothing of the well-studied self-adjoint 
case flLSjl . 

1.3. Layout of the paper: duality of global and blow-up asymptotics. In Section 
[21 we describe some properties of the fundamental solution of (II. ip given by 
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(1.10) bix,t) = t~2^Fiy), y=^ =^ BF = in 

Next sections are mainly devoted to the following two asymptotic problems for (II. ip : 

Application I: Global asymptotics as t -)■ +oo. Sections [3] and [51 The asymp- 
totic behaviour as t — )■ +oo of solutions is governed by the eigenfunctions of (II. 6p : 

(1.11) B = -i(-A)- + ^yV + ^/, witha(B) = {A, = -M, |;j|>0}, 

which demands a proper definition of its domain by a spectral decomposition (Section 
[3l) and by a traditional spectral theory involving careful using poles of the resolvent, 
non-classic "radiation conditions" at infinity (see an alternative approach in [26]), etc. 
More precisely, we establish that the discrete spectrum and the eigenfunction set for the 
operator (II. lip describe all the possible asymptotics as t — )■ -|-oo of solutions of (II. ID for 
any data Uq G L^*(]R^). The exponential weights (recall the change p \-^ p* relative to 
(II. 8p done by some clear and natural reasons) 

(1.12) p*(y)=e\y\\ with a = ^, and p(y) = -^ = 6"^", 

are properly introduced in Section [31 It is curious that, even in the classic case m = 1, 
i.e., for (II. 2p . we were not able to find any essential traces of such a full refined scattering 
theory (except some particular results often admitting not-that-clear interpretation) and 
corresponding spectral properties in the vast existing scattering literature!. 

The classic real analogy and forerunner of (II. lip is the self-adjoint operator for m = 1 
(we apologize for the necessary change of the weights here, p ^-^ p*; cf. (I1.12p ) 

(1.13) B = A + i^y-V + f / = iV-(pV) + f / in ^^(M^), with p = e^. 



The authors do not stiU beUeve that optimal and sharp large-time (t — >■ oo) asymptotic theory for the 
classic LSE (|1.2I) in M^ has not been developed in full details (some particular results have been indeed 
known seems) since this PDE burst into quantum mechanics and mathematical physics in the 1920s. As 
was mentioned, a similar (and even stronger) asymptotic theory for the ID heat equation ut — Au (cf. 
(|1.1|) for 777, = 1) is known since Sturm's analysis j85' of formation of multiple zeros of solutions (we will 
do this for ()1.1|) in Section[6]) obtained in 1836! The authors will be very pleased to get rid of such a quite 
surprising delusion, but also would be naturally satisfied to know that this, though again unbelievably, 
is done in full details for the first time in the present paper. 
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As we have pointed out already, its real discrete spectrum cr(B) = { — 2, ^ > 0} and 
eigenf unctions as Hermite polynomials, multiplied by the Gaussian, are known from, at 
least, the 1830s, and are associated with the names of Sturm and Hermite; see |23l § 1.2] 
for more history and original Sturm's calculations, and pi p. 48] for a fuller account of 
applications of these separable polynomials in self-adjoint linear operator theory. 

Thus, we are obliged here to develop Hermit ian-like spectral theory for the 2mth-order 
rescaled Schrodinger operator (11. lip , and this is an unavoidable task if we want to reach 
an optimal classification of large time behaviour for the non-stationary LSE (11. ip . 

Application II: Blow-up asymptotics as t -)■ 0~, Sections m El Appendices B 
AND D. Alternatively, for data Uq G L^*(M^), using blow-up scaling at a finite point as 
a; — )■ and t — )■ T^ = 0^, we show that this behaviour of solutions and local structure of 
their nodal sets are described by the eigenfunctions of the linear operator (II. 4p . which is 
"adjoint" to B in a sense (but not in the standard dual L^-metric; an indefinite metric 
should be involved, which will be carefully explained). The discrete real spectrum (I1.7P 
remains the same as in (II. lip . 

A key point is that the eigenfunctions of B* are generalized Hermite polynomials 
{ip*^{y)}, so that the nodal sets of solutions of (II. ip are locally governed by zero sur- 
faces generated by these polynomials only (there exists a countable, complete, and closed 
set of those). Application III: This allows us to state in Section [6] a sharp uniqueness 
continuation theorem for (II. ip : if lrsiu{x,t) (or Iieu{x,t)) has, roughly speaking, 

a rescaled nodal set component not decomposable into a finite combination 
from polynomial surfaces {liatpl{y) = 0, |/3| > 0}, then u(x,t) = 0. 

Some further applications of these spectral results are also discussed. 

In Sections \5\ and El possible applications of Hermitian spectral theory are studied for 
more general linear PDEs and for the 2r?T,th-order nonlinear Schrodinger equation (the 
NLSE-2m) 

(1.15) Ut = -i{-A)'^u±i\u\P-\ in M^ x M+, where p>l 

and the sign "-(-" corresponds to the focusing (blow-up) model, while "— " gives a defo- 
cusing one. See [HI |66l |671 |77l [811 193] as a guide concerning the modern research of both 
semilinear PDEs (ll.lSp . 

Application IV: In Appendix B, we show how to apply the spectral results to the clas- 
sic problem of the regularity of a boundary characteristic point for the linear Schrodinger 
equation. It seems that, in the present setting of rather arbitrary "backward paraboloids" 
at characteristic points, such issues were not addressed in the existing literature. 

(ii) Application V: Quasilinear Schrodinger equation, Appendix C. This is a 
most "risky" application of our refined scattering spectral theory for the spectral pair 
{B,B*}, so we also put this into an appendixlj- Thus, exhibiting a certain necessary 



This and some other non-entirely-rigorous parts of our apphcations well resonate with two Tao's 
comments in his "What is good mathematics?": (i) "... mathematical rigour, while highly important, 
is only one component of what determines a quality piece of mathematics" [571 P- 624], and (ii) "... we 
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bravery, we develop some basics of a "nonlinear eigenfunction theory" for a quasilinear 
Schrodinger equation (the QLSE) of the form 

(1.16) Mt = -i(-A)"'(|M|"ii) in M^ x M+, where n>0 

is a fixed parameter. The applications and some history/references concerning such rather 
unusual quasilinear PDEs are explained therein. Here, we intend to reconstruct a proper 
connection between linear and "nonlinear" spectral theory by performing a continuity 
homotopic path n — )■ O"*", which establishes a link between "nonlinear eigenvalue problems" 
for (11.161) and the linear one developed for (II. ip . As a result, we predict existence of a 
countable family of the so-called n-branches of solutions, which are originated at n = 
from eigenfunctions on the corresponding eigenspaces for the linear spectral pair {B, B*}. 

We strongly believe that the results of Hermitian spectral theory developed can be 
useful for attacking a number of open problems concerning blow-up behaviour for (ll.lSp . 
(11.161) . and others. We plan to explain this in a forthcoming paper. Meantime, we just 
comment on that standard blow-up rescaling {q.v. (14.61) ) leads to the adjoint operator 
(II. 4p as the linearization, so the generalized Hermite polynomial eigenfunctions of B* can 
be key for understanding this intriguing rescaled blow-up dynamics. 

2. Fundamental solution and the convolution 

2.1. Fundamental solution and its first properties. In constructing fundamental 
solutions and corresponding convolutions for (II. ip . one can use the fact that formally 
changing the independent time variable it \-^ t yields the standard poly-harmonic PDE: 

(2.1) it^t =^ ut = -i-A)"'u in M^ x M+. 

This creates an artificial complex (imaginary) time axis, and nevertheless will help us 
to restore various spectral properties and other functional details related to (II. ID . Of 
course, unlike the real parabolic case [17], the change in (12.11) implies a well known highly 
oscillatory properties of the fundamental and other solutions of the Schrodinger equation 
(II. ip . that are not available for its real parabolic counterpart in (12.10 . 

Thus, by classic PDE theory, given proper initial data Uq{x), the unique solution of the 
Cauchy problem for (12. ip is given by 

(2.2) u{t) = b{t) * Mo for t > 0, 

where b{x,t) is the fundamental solution (ll.lOp of the operator in (12. ip . Substituting 
b{x,t) into (II. ip . one obtains the rescaled kernel F{y) as a unique solution of a complex 



should also be aware of any possible larger context that one's results could be placed in, as this may 
lead to the greatest long-term benefit for the result, for the field, and for mathematics as a whole." [87l 
p. 633, Tao's last sentence therein]. [In author's opinion, Tao's (ii) statement could definitely serve as 
a most impressive characterization of a "good mathematics" among dozens of others presented in his 
paper.] The latter one (ii) somehow helps us to justify including some "risky applications" , which we are 
not aware how to prove, and even are not sure whether such results can be proved in any remote future; 
though the authors believe that this stuff should be revealed for the Readers, which might be interested 
in nowadays or will be later on. 
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linear ordinary differential equation (ODE) which is the radial restriction of a linear PDE 
system, 

(2.3) BF = -i{-^yTF+^yWyF+^F = Q in M^, 

so there occurs the linear operator B given in (II. lip . In addition, the kernel F is defined 
in such a way that, in the sense of distributions (or other suitable L^-type metrics): 

(2.4) 6(t)*no^Mo, t^0+, or j F{y)x{t^y)dy ^ xi'^) Vx G Co~(M^), 

justifying initial data. These define the unique rescaled kernel F. 

On the other hand, using the Fourier transform yields the following equivalent repre- 
sentation of F\ 



(2.5) HK-.m^) = e-^l^l ""* =^ HF{-)){y) = e-'l^l ""• 
For m = 1, this gives the "Gaussian" exponential profile 

(2.6) F(?/)=(^-i^e^ (m=l). 

It follows that F{y) is highly oscillatory a.s y ^ oo. In particular, for (12. 6p . we have 

(2.7) F(y) = (^-i^[cos(M!)+isin(M!)], |F| ^ const., F ^ D>{R^), p>l. 

For arbitrary m > 2, the asymptotic behaviour of F{y) for \y\ ^ 1 is covered by the 
classic WKBJ asymptotics. Namely, fixing in (12. 3p two main leading terms for the radial 
kernel F = F{z), for z = \y\ —> oo, 

(2.8) -i(-l)'^F(2"^) + 2^F'z+... = 

yields, in the first approximation, the following exponential asymptotic behaviour: 

(2.9) F{y)^e^\y\" =^ « = ^1^ and (aa)^™-! = ti£^. 

Hence, there exist 2m — 1 different complex solutions {ofc} belonging to a circle in C: 

(2.10) Ittfcl = Zm = ^ (2m)"2^^^^ < 1 for m > 1. 

Obviously, we are interested in those roots Ofc, for which Reafc < 0. Otherwise these will 
be exponentially growing oscillatory solutions that will be "too much" non-integrable. It 
is clear that there exists the purely imaginary root with the main asymptotic oscillatory 
behaviour at infinity: 

(2.11) ao = Zmi =^ F{y) r^ cosi^Zmlyl") + isin{zm\y\°') as y ^ oo. 
On the other hand, the ODE (12.31) admits solutions with a power decay: 

(2.12) J-F'z + ^F + ... = 0, 2=|y|»l =^ ^{y) ^ ^ as y ^ ^. 

Since F{y) is then are non-oscillatory (for being used as in (12. 4p ) and, in addition, are "too 
much" non-integrable as ?/ — )■ oo, such asymptotics are not acceptable for the fundamental 
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kernel F{y). We then arrive at the following simple, but interesting and, in fact, a key 
property of F{y), which will affect our analysis (especially, in the "nonlinear" cases): 

(2.13) all asymptotic components of the rescaled kernel F{y) as y — )■ oo are oscillatory. 

2.2. Convolution: a unitary group. Thus, the unique weak solution u{x, t) of the 
Cauchy problem (II. ip for any data Uq G C is given by the Poisson-type integral for t eM: 

(2.14) u{x,t) = b{t)*uo = e-'^-^^"''uo = t-^ J F{{x - z)t-^)uo{z) dz, 

where {e~'^~^^™*}ig]R is the corresponding unitary group. In what follows, for some con- 
venience, we take t > only (that implies no trouble in defining the flow in (I2.14p ). so 
actually we deal with the semigroup {e~^^~^^"'^}t>o- 

3. Discrete real spectrum and eigenfunctions of B 

This section is devoted to some preliminary analysis of spectral properties of the key 
pair of linear rescaled operators {B,B*} that appear after long-time {t — )■ +oo) and 
short-time (t — )■ T~) respectively rescaling of the LSE-2m (II. ip . In fact, this explains in 
a reasonably brief manner several necessary key properties of the pair to be used later on. 
However, we must admit that some of the issues will require a hard work to justify by 
classic theory. This will take a full separate paper |2S] . Nevertheless, we hope that listing 
a full collection of some involved spectral properties will be convenient for at least some 
of the Readers, who are interested in general understanding of how this approach works 
and who do not require full mathematical details. 

3.1. First step to the domain of B in a weighted L^-space. We now study spectral 
properties of the first appeared linear operator B given in (12. 3p in the space Lp(M^) with 
the exponential weight: 

(3.1) p{y) = e-\y\" > in M^, where a = ^^. 

By (■,■), we denote the standard L^-product: 



(3.2) {v,w) = J v{y)w{y)dy. 



As customary, if^™'(R^) denotes a Hilbert space of functions with the inner product 



p 

2m 



(3.3) {v, w)p = J p{y) E D^v{y) D^w{y) dy, 

RiV fc=0 

where D^v stands for the vector {D^f , |/3| = A;}, and the norm 

2m 

(3.4) \\v\\l= j p{y)Y.\D'v{y)\''dy. 

RiV fc=0 

T2m/ia,N\ r- T2fTn>N\ -^ T2fjn,N 



Obviously, if;'"(M^) C L;(M^) D L\R^). Introducing the weighted Sobolev space 
H'^"^{M.'^) is a first step to better understanding the necessary and natural domain of 



B, as stated in the proposition below. However, the space L^(]R^) with the exponen- 
tially decaying weight (13. ip is evidently too wide, so we cannot expect any good spectral 
properties therein. Nevertheless, we now prove the following: 

Proposition 3.1. B is a bounded linear operator from i7^™'(R^) to L^(M^). 
Proof. It follows from (Q that Bv e ^^(M^), if 

(3.5) J p{y)\yVv\'^dy<C\\vfp for any v G //^m^^^V)^ 

where C > is a constant. Let h G C°°(R^) be a function such that h{y) = 1 for \y\ > 2 
and h{y) = for \y\ < 1. Since the inequality 

J p{y)\yV[{l-h{y))v]\'dy<C,\\v\\l 

is obvious, it suffices to show that 

J piy)\yVihv)\'dy<C\\hv\\% 

R^ 

i.e., proving (13. 5p . we can suppose that v G C^(M^) vanishes for all \y\ < 1. 

Let (r, 6*1, ..., ^Af_i) be the spherical coordinates in M^. Since \y ■ Vv| < r\vr\, it suffices 
to verify that 

oo oo 

(3.6) / r^+Vi(r)| V^" dr < C^J r^- Vf™~'V)l'e-^" dr, 



if the left-hand side in bounded, and apply this estimate with wi = Vr- 

Let (3' = A^ — lorg = A^ — 3, 7 = a — 1 = ^m-i ^^^ ^ ^ (0' 1) ■••) 2m — 2}. Then using 
the inequality 

oo 

J ri+^k^\w'{r) + r^w{r)\'^e-''" dr > 0, 



integrating by parts again implying that the right-hand side converges, we obtain that 

oo oo 

Jri+^k^\w'{r)\^e-'^"dr > 7 / r'?+2'=^+2^|w(r)|2e-^" dr. 


Simple iteration implies fl3.6p with C2 = ^i"^'", completing the proof. D 

The result also follows from a general estimate in [391 Lemma 2.1], which goes back to 
the Hardy classical inequality [37]. In a similar (or obvious in (ii)) manner, introducing 
the "adjoint" spaces with the reciprocal weight LL(M^), with the weight fl3.7p . 

(3.7) p*(y) = ^ = eH^ 

we have the following: 

Corollary 3.1. B is bounded as an operator 
(i) B : //2m(^iV) ^ L2,(M^) and 

(ii) B: Hf^{R^)^L%R^). 
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Remark for m = 1. As customary, in the second-order case 171 = 1, there appear 
some extra possibihties and "symmetries". Namely, then B admits a formal symmetric 
representation 

(3.8) B = i i V ■ (kV) + f /, with the "weight" K{y) = e"' ^^ , 

in the weighted space L^ with the complex weight k ^ and hence with an "indefinite 
metric"; cf. Azizov-Iokhvidov |1], which we will need to refer to later on a few times 
at least. We do not know any reasonable application of the complex symmetric form 
in (13. 8p . For instance, as usual, the symmetry (13. 8 p implies the formal orthogonality of 
eigenfunctions: 

(^/3, ^7)« = J ni/jpip-y = for /3 7^ 7, 
but it is not that easy to find suitable applications of this in view of the indefinite metric 
involved. Anyway, by no means, we are going to rely on this kind of a pseudo-symmetric 
representation of the operator for m = 1, especially, since for m >2, this illusive complex 
symmetry of B disappears without a trace. 

3.2. Group with the infinitesimal generator B. Before introducing detailed spectral 
properties of B, we present a simple derivation of its group for proper weak solutions to 
be heavily used in what follows. 

Thus, the rescaled solution of (II. ip defined as 

(3.9) w{y,T) = t'^u[yt2^,t), where r = lntGM (t > 0), 
satisfies the necessary rescaled equation 

(3.10) Wr = B-u; (the operator B is as in (12. 3p ). 

Then w{y, r) solves the CP for (13101) in R^ x M+ with data at r = (i.e., at t = 1) 

(3.11) Wo{y) = u{y, 1) = b{y -■,!)* %(■)• 

Rescaling convolution (I2.14p yields the following explicit representation of the group with 
the infinitesimal generator B: 

(3.12) w{y,T) = e^^u{y,l)= J F{y- ze-^^)uo{z)dz for t eR. 

Performing another rescaling 

(3.13) w{y,T) = {l+t)^u{y{l + t)^,t), where r = ln(l + t) : M+ ^ R+, 

we obtain the solution w{y, r) of the Cauchy problem for equation (I3.10p with initial data 
Wo{y) = Uo{y). Rescaling (12.140 . we deduce a more comphcated, but standard (without 
the relation (13.1 ip ) representation of the semigroup for r > 

(3.14) w{y,T) = e^^uo = (1 - e"^)"^ / F{{y - ze-^^){l - e-^)-^)uo{z) dz. 

By the Holder inequality (see e.g., (I3.24p below), it is easy to see that 

(3.15) w{-,t) e Ll{R^) for all r>0 (mq G LJ,(M^)). 
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3.3. Spectral decomposition of B: a first step to discrete spectrum and eigen- 
functions via converging expansion of convolution. We now in a position to make 
more clear a proper definition of the necessary operator B to be exploited later on. We 
then confirm the actual existence and the so-called evolution completeness/closure of some 
eigenf unctions for initial data uq G L'j^,(M'^). It is worth to stress out how simple such 
derivations are here, unlike a more standard spectral and analytic continuation techniques 
applied in [26] to justify necessary hard properties of resolvent poles and related issues. 

It can be derived from the ODE (12. 3p for the rescaled kernel F{\y\) and also from (12.51) 
that the higher-order derivatives of F can be estimated as follows: 

(3.16) \D^F{y)\<C{l + \y\) 



(-1)1/31 in M^ 



Actually, according to (12.91) . (I2.10p . a sharper estimate includes the factor |a2;m,|^ i-e., 

(3.17) la-^ml' = (2m)~2^^^^ -s. o as / -)■ oo, 

so this improves convergence of the series to appear later on. 

Consider Taylor's power series of the analytic kernel F{-) on compact subsets y ^ {\y\ < 
L}, with an L ^ 1, 

(3.18) F{y - ze-^) = Ee"!? ^-^D^F{y)z^ ^ E^~'-^^.My>', 

(/3) W 

where z^ = z^^...z^ and ipis are in fact normalized eigenf unctions of B; see below. This 
series converges uniformly on compact subsets in 2; G M^. Indeed, for |/3| = / ^ 1, we 
have the following approximate estimate of the expansion coefficients: 

(3.19) E ^-j^D^Fiy)zf\..z^/ < C}, (1 + \y\)(-~'Hl + \z\y, 

where we have used a rough bound by Stirling's formula (as usual, we often omit the 
lower-order multipliers ~ C' in (I3.19P ): 

/3!>[(^)!]^~iV'/!. 
Finally, we arrive at the following representation of the solution: 

(3.20) w{y,T) = Ee-^M^(«o)V'/3(2/), 

where -^/^ = ~2m ^^^ ''Ppiv) ^^^ ^^^ eigenvalues and eigenf unctions of B and 

(3.21) Mp{u,) = ^ j zP,\.J^-u,{z)dz 

are the corresponding moments of the initial datum wq (recall the relation (13. lip between 
wq and Uq). We will show that, in terms of the dual inner product (■, ■) in L'^{EJ^)., 

(3.22) M^(mo) = K,V';), 
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where {i^l} are polynomial eigenfunctions of the adjoint operator B* to be described in 
greater detail in Section HI It is not difficult to check that ( I3.20p uniformly converges on 
any compact subset y E {\y\ < L}, since for / = |/3| ^ 1, 

(3.23) \M^iuo)My)\<C}^,L(--'^' Jil + \z\y\uoiz)\dz, 

where we estimate the last integral by using the radial variable and the Holder inequality, 

/ (1 + \z\y\uo\ dz= f{i + \z\y-^ ^\uo\ dz 

(3.24) ^"^ "^ 

< \//iRiv(l + kl)^'irdz Jf^^p*\uo\^dz. 



p 

The last integral is bounded since uq G L^* by the assumptions, while the ffist integral 
for very large / can be roughly estimates by Stirling's formula as follows: 

(3.25) v//m.(1 + \z\rj, dz ~ ^/- r^^-n-^'e-- dr ~ (^) ^ ~ {^f^li. 

Since a > 1, the right-hand side is essentially smaller than /! ~ /' in the denominator in 
(I3.23p . Hence, this series converges uniformly and to an analytic solution as expected, 
and not surprisingly for the LSE. 

Recall that (I3.14p gave the actual corresponding group {T{r) = e^'^j^gR. Here, we can 
apply the same expansion analysis as above, which directly determines the eigenfunctions 
of B*. First, by Taylor's expansion in the F-term, we obtain 

(3.26) w{y, t) = J2 ^^^D^'F{y{l - e-^)'^) (e^ - l)-^^ / z^'woiz) dz. 
Second, applying Taylor expansions for the /-terms in (I3.26p . 



(3.27) F{yil - e-)-!^) = ^h hD-^FiO)y^ (1 - e^)"^ 



N . „ kT 



and for functions (1 — e '^) 2m and (1 — e '^) 2m in terms of e 2m , A; = 0, 1, ..., we arrive 
at a similar representation of the semigroup 

(3.28) w{y, t) = e^-wo = E(/3) e''-^ Mp{wo)My)^ ^ > 0, 

where the expansion coefficients Mi3{wo) are dual products (tocV'^) with the polynomial 
eigenfunctions f/^S of the adjoint operator B*. Expansion (I3.28P determines the adjoint 
eigenfunctions {■0^} that are "orthogonal" to {"0^} in L^(]R^) (all products and metrics, 
with the convention (14. 4p . to be introduced later on). In Section H] devoted to the adjoint 
operator B*, we perform a simpler derivation of explicit formulas for polynomials {i^t}- 

Thus, using the expansion (13.280 . we are now in a position to present the ffist proper 
definition (via its spectral decomposition) of our bounded linear operator in i7^™(M^), 
which we denote by B: 

(3.29) Mwo = £e*^wol^^ = j:i^)\0Mp{wo)My) for woeH^^n^^). 
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Actually, IB can be considered as a restriction of a more general operator B defines in a 
so-called wider space of closures. Since this extension is not essential for our main appli- 
cations and B is sufficient, we postpone this rather technical procedure until Appendix D 
at the paper end. 

Comment on "extended eigenfunctions" . With the definition of our operator 

(3.30) B : if^l'"(M^) ^ ^^(M^) 

by its expansion f l3.29p . we face the following technical difficulty. Namely, the eigen- 
functions {"0^}, which actually generate such a B, do not belong to its domain in f l3.30p . 
and we refer to them as to extended eigenfunctions (i.e., as we will show in Appendix 
D, belonging to an extended space). However, for simplicity, we continue to call them 
simply eigenfunctions, bearing in mind that, in Appendix D, we are going to construct an 
extended space of closures, where all ifj^ belong to. We then restore the original operator 
B rather than its restriction B, though, as we mentioned, for our main PDE applications, 
present "spectral theory" of the restriction B is more than sufficient. 

We now summarize the above results concerning the introduced operator B. Recall 
that, originally, B is defined by the rescaled convolution (I3.14p . with the corresponding 
space and the domain. For further applications, we restricted B and defined it via the 
spectral decomposition such as (I3.29p . which demanded special topologies for a proper 
convergence. 

Proposition 3.2. The semigroup expansion series ( \3.2^ for wq G LL(]R^), and the 
eigenfunction expansion o/B in i\3.29[) for wq € HpT'lM.^) converge: 

(i) uniformly on compact subsets in y, and 

(ii) in the mean in L^(]R^). 

Proof, (i) has been proved, (ii) The convergence in the mean in L^(]R^) is not straightfor- 
ward for such a bad oscillatory and growing basis functions in $. Therefore, estimating 
the terms in 03.2Up . we need to include both decaying multipliers in the right-hand side 
of (I3.25p . so that the following estimate of a typical integral is key by using Stirling's 
asymptotics of the Gamma function: 
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(3.31) "° "° 

Since a G (1,2], this exponentially (but not superexponentially as before for the coeffi- 
cients of (I3.20p on compact subsets in y) decaying coefficients guarantee the convergence 
of the series in L^(R^). The factor in (I3.17P is then not necessary for convergence. D 

3.4. Discrete spectrum. Thus, the series (I3.29P is a spectral decomposition (actually, the 
eigenfunction expansion for a discrete spectrum) of our linear non self-adjoint operator, 
which we have denoted by B. Hence, we treat it as a linear bounded operator given in 
(I3.30p . when the eigenfunction expansion (I3.28P is defined in two topologies: local uniform 
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convergence, and strong convergence in the mean in L'l{W'^). In this case, the expansion 
in L^»(M-'^), with the same convention on topologies (see (14 ■4p explaining why if^l occur 
in the last products), 

(3.32) weL%{R^) =^ w = Zcf,i^p, cp = {w,rp), 

is naturally treated as the eigenfunction series representation of the embedding operator 
L2,(M^) -^ L^iR^) in terms of the eigenfunctions $ of B. 

For convenience and more systematic understanding the new class of linear operators 
introduced, we clearly state necessary spectral properties of the operator B. According to 
[26] . the present operator B, which was not still properly defined and has been uniquely 
characterized by its "spectral eigenfunction decomposition" (13.291) only, is not allowed to 
have the positive part of the point spectrum with polynomial eigenfunctions, as shown in 
[26], where this is done by a direct introducing "radiation conditions" posed at infinity. 

Thus, according to our operator definition (I3.29p . we ascribe to B a countable set of 
numbers, which, for convenience, we continue to call its discrete spectrum: 



a{M) = {X^ = -^, 1/31 = 0,1,2,...}, 



(3.33) 

where eigenvalues A^ have finite multiplicity with eigenfunctionci 

(3.34) 



My) = ^D'm - ^-^a-r-i^Y^ny). 



The existence of such eigenvalues and eigenfunctions is dictated by (13. 29 p . It is worth 
mentioning that the same follows by applying D^^ to the elliptic equation (12. 3p (here B 
stands for its differential expression, so we are not obliged to use B): for any f3, 

(3.35) L"^BF = BD^F + U D^F = =^ {D^F, A;, = -g} is a pair for B. 
Let us fix some other properties of such extended eigenfunctions: 



Lemma 3.1. (i) The subset of eigenfunctions $ = {ip/s} is complete in L'^iM.'^), and 

(ii) $ is L'^p, -evolutionary closed in the sense that the eigenfunction expansion {\3.28^ . 
which converges in the means and uniformly on compact subsets, presents the rescaled 
solution of the LSE (IJ. jp for any data uq G L'^*{EJ^); and 

Remark: towards more general integral evolution equations and rescaled op- 
erators B. As we have mentioned, all the above results can be justified by classic spectral 
methods associated with the given rescaled differential operators, |26j. 

However, it is worth mentioning now that all the (i)-(iii) remain valid for more general 
class of integral evolution (pseudo-differential) equations (I3.12p . where 



(3.36) F{y) is an arbitrary sufficiently good analytic kernel. 



In other words, (13.120 then do not correspond to any linear PDE. Furthermore, the 
expansions (I3.18P and (I3.20p can be also prescribed. Overall, this gives the spectral 



Actually, extended eigenfunctions, since "0/3 ^ i/pI"(]R^); a standard meaning eigenfunctions of B to 
be restored in Appendix D by introducing an extended space of closures of finite eigenfunction expansions. 
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results similar to those in (i)-(iv), where extra efforts to justify the functional topology 
required are necessary. 

In this connection, it is key to emphasize that, in the present most general situation, 
no extra powerful tools of spectral theory, developed in [26] for the present Schrodinger 
operators, will be at hand. Then, we will be inevitably attached to a different functional 
framework, with no visually available operator B (and the "adjoint" one B*). 



Returning to Lemma 13. ![ note that, as an important characterization of the eigen- 
function set $, due to (12.131) . all of the eigenfunctions il'/3{y) satisfy the property (I2.13p . 
i.e., these do not have any "non-fast-oscillatory" asymptotic component as y — )■ oo (since 
(12.121) has been excluded from the kernel F[y)). Observe also from (13.341) that, with a 
proper definition of such extended (generalized) linear Junctionals, to be done in Appendix 
D, the following can be interpreted as being correct: 

(3.37) (V^o, ^o)* = / My) dy= J F{y) Ay = 1, since Vo = F, ^o = 1- 

Recall that, in the usual sense, such oscillatory integrals are not properly defined. On the 
other hand, it also follows from (I3.34p that, again in a proper extended linear functional 
sense (see Appendix D), 

(3.38) (V'/3,V'S)* = / i^M^y = for any \(3\ > 1. 

This is easier to believe in view of the integration by parts, though the oscillatory inte- 
grals are not well defined as well. In fact, these equalities express the orthogonality of 
any i/jp to the first adjoint eigenfunction ip^ = 1 via the dual inner product, defined as 
extensions of linear functionals prescribed in (I3.32p : see again Appendix D. The adjoint 
eigenfunctions are polynomials which form a complete subset in L^(]R^) with the same 
decaying exponential weight (13.11) : see Section HI 

In the second-order case m = 1, using the rescaled kernel (12. 6p in (13.341) gives the cor- 
responding generalized Hermite polynomials Hp{y) (given up to normalization constants) 
via the generating formula: 

(3.39) My) = jdymD^^''-^ = Hp{y)e''4. 

Note that, in this case, the generalized Hermite polynomials are obtained from the classic 
ones [U p. 48] by the change 

Hp{y) = Hf'-^-{\y) {m = l). 

Proof of Lemma \3.1\ 

(i) Completeness in LFp-space. In order to prove completeness in the metric of L^, as in 

[T71 § 2], we suppose that there exists some function G (say, G G L^), which is orthogonal 
relative to the inner product in L^ to all eigenfunctions, i.e., 

/ p{y)D''F{y)G{y) dy = for all a. 
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Since F is analytic, it implies that 

J p{y)F{y-x)G{y)dy = for all x G M^. 

Consider the Cauchy problem for (11.11) with initial data 

uo{x) = p{x)G{x) in R^ (note that pG G LJ.(M^)). 

One can see from the Poisson-type integral fl2.14p that the solution exists for all t G (0, 1]. 
Then u{x, t) is analytic in x. We have 

u{x, 1) = / F{x - y)G{y)p{y) dy. 

Therefore, u{x, 1) = 0. It follows from the standard uniqueness theorem (see [3S] as a 
guide) that u{x, 0) = 0, and G = 0. 

(ii) Evolution closure in L^, has been already proved while studying the convergence of 
the series (K^ and flX^ . 

This completes the proof of Lemma 13.11 D 
Remark: positive part of the "point spectrum". Formally, in addition to ([ 



there exists another positive part of the "real spectrum" of the differential form B (12. 3p : 

= {^/^ = ^' 1/31 = 0,1,2,...}, where 



(3.40) 

i^piv) = 'JW.y^ "*" ■•■ ^^^ |/3|th-degree polynomials. 

We will show how to construct such polynomial eigenfunctions in Section HJ where these 
are actual eigenfunctions of the adjoint operator B*. Though such "pseudo-Hermite" 
polynomials $+ = {-i/'t G L^} actually exist, they are not eigenfunctions since are not 
available in the given expansions (I3.28P as the definition of the operator. In other words, 
these do not satisfy the necessary radiation conditions at infinity, which turn out to be 
rather non-standard and unusual, [261. 



4. Spectrum and polynomial eigenfunctions of the adjoint operator B* 

4.1. Indefinite metric and domain of the bounded operator B*. Using the results 
obtained above for differential expression B and its proper restriction B in (I3.29p . we now 
describe in detail the eigenfunctions of a restriction of the "adjoint" operator (ll.4p . which 
is still a differential expression and we use the notation B and not B*. B* will be obtained 
via blow-up rescaling (14. 6p . Note that B* is not adjoint in the standard dual metric (13. 2 p 
of L^, since, as we have seen, 

(4.1) B* = (B)*^, = -B + ^/ = B*. 

Curiously, B* is then the standard adjoint to B (as a differential form, in C^(M^)) in 
the indefinite complex metric (cf. (13. 2p ) 

(4.2) («, w), = J v{y)w{y) dy =^ (B): = B*, 
where the complex conjugation in the second multiplier is not assumed. 
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In fact, being treated in the metric of L^, without a weight as it used to be, such an 
indefinite metric is not that challenging and even can be partially get rid of for the practical 
use of eigenfunction expansions; see below. However, some comments are necessary. 
Firstly, the set of real L^-functions 

E+ = {veL'^: lmv = 0} 

is a positive lineal (a linear manifold in the field of real numbers)) of the metric, i.e., 

{v, f )* > for V e E+ d L^ , V ^ 0. 

The purely imaginary functions, 

E_ = {veL^: Rev = 0}, 

define the corresponding negative lineal. Therefore, L^ with this metric is decomposable: 

v = v+ + v. = '-^^^^p^ + '-^^^^^, where v± e E± ^ L^ = E+ ® E_. 

This defines the corresponding positive majorizing metric as follows: 

\{v,v)^\ < [v,v]^ = {v+,v+)^ - {v-,v^)^, 

etc. It should be noted that such a case of the decomposable space with an indefinite 
metric having a simple majorizing one is treated as rather straightforward; see Azizov- 
lokhvidov Pj for linear operators theory in spaces with indefinite metricqj. 

Though, as we have mentioned, linear operator theory in spaces with indefinite metrics 
exists for more than half a century, we do not think that the complex indefinite metric 
in f l4.2p creating the necessary pair {B,B*} (the operator and its adjoint) can play any 
role in what follows. On the other hand, as customary, using the metric f l4.2p is not that 
suspicious, since it is necessary only for calculating the expansion coefficients according 
to the standard rule: 

(4.3) v = ^Ci3'iIJ(3 =^ C,3 = {V^lpl)*, 

while all convergence calculus can be continued to be performed in standard metrics. 
However, to avoid possible future accusations of using non-approved indefinite metrics, 
we are now back to standard scalar products by noting the following. Since B* is shown 
to have a real point spectrum only, using the standard L^-metric instead of (14. 2 p will only 
mean replacing the eigenfunctions as follows: 

(4.4) ^^^v^;, 

and we are assuming using this convention any time when necessary and convenient. 



Basic results of linear operator theory in spaces with indefinite metrics can be found in Azizov- 
lokhvidov's monograph jj. It was in 1944, when L.S. Pontryagin published the pioneering paper "Her- 
mitian operators in spaces with indefinite metric" j78j . A new area of operator theory had been formed 
from Pontryagin's studies, which, during the time of the WWII, were originated and associated with 
some missile-type military research (a comment by Yu.S. Ledyaev). This work set by Pontryagin was 
continued from 1948 and in the 1950s by M.G. Krein [521 [SD], I-S. lokhvidov gD], and others. 
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In the second-order case m = 1, (11.41) has a formal complex symmetric representation 

(4.5) B* =i^V-{K*V), K*(|/) = e'T^, 

though we do not use this. Similar to B, we do not know any advantages, which this 
symmetry in such an indefinite metric can provide. However, for m > 2, any formal 
additional symmetry is not available. 

For any m > 1, we again consider B* in the weighted space L^(M^) with the same 
exponentially decaying weight (13. ip . and ascribe to B* the domain Hp"^(R^), which is 
dense in L^(M^). As in the previous section, B* : if^''"(M^) — > L'j^(R^) is shown to be a 
bounded linear operator. 

4.2. Semigroup with infinitesimal generator B*. In order to construct the semigroup 
with the infinitesimal generator B*, we use the rescaled variables corresponding to blow-up 
as t — >• 1~, 

(4.6) u{x,t)=wiy,T), y = ^^-^, T = -ln(l-t):(0,l)^M+. 
Then w solves the problem 

(4.7) Wr = B*w forr>0, w{0) = Uq E Ll,{R^). 

Rescaling (I2.14p . we obtain the following explicit representation of the semigroup: 

(4.8) w{y,T)=e^*^uo = {l-e-^)-^ J F{{ye-^^ - z){l - e-^)-^)uo{z) dz. 

4.3. Spectral decomposition and definition of B*: using explicit representation 
of the semigroup. Similar to B in Section |3l the original rescaled "adjoint" operator 
B* is defined by the convolution (14. 8p . For the purpose of applications, we will need its 
restriction defined in terms of its spectral decomposition obtained via the semigroup rep- 
resentation (14. Sp . Comparing semigroups (14. 8 p and (I3.14p . we see that the only difference 
is in the argument of the rescaled kernel F{-). Therefore, instead of (I3.27p . we have to 
use the following expansion: 



F((y e-^^ -.)(!- e-r^)=Ew^^(l/e-^^- ^)'(l-e-r- 



(4.9) 

where (?/e-2^^ - 2)^ = ^(o<5<^) C^e'^^ y^-^ {-zY . 

Then, using both ( KM in (SSI) yields 

(4X0) (^>0) (7) (0<5<7) 

xi-iPy^-^..{l±^)j^Dl^[J{l,D^FiC)z')uoiz)dz]\^^^. 
This is the expansion over the point spectrum of B*, 

(4.11) w{y,T) = e^'^wo = j:e-^^M*{u,)rp{y), 

W 
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where ip*a{y) are finite polynomial eigenf unctions (see their direct derivations below) and 
the expansion coefficients are 

M*(uo) = (mo,'0/3)- 

Convergence of the series fl4.1ip is studied as in Proposition 13.21 Similar to f l3.29p . the 
group representation (14.111) defines the necessary operator B* satisfying (I3.30p as 

(4.12) Wwo = i^e»*-w;o|^^o = ^ A/j M^(ti;o)^^(2/), Wo e i/2-(R^). 

(/3) 



Similarly to the case of the operator B at the beginning of the proof of Lemma 13. H in 
view of standard regularity properties of linear parabolic flows such as (14. 7p . the semigroup 
expansion (14. lip reveals some key auxiliary spectral properties of B*: 

(i) the point spectrum is a{M*) = {A/3 = — ^}, with any A^ having finite multiplicity; 

(ii) by the definition, there is no continuous spectrum; and 

(iii) polynomial eigenfunctionq^ {'^^(y)} are closed in L^, etc. 

In addition, we have to observe that, unlike (I2.13P for B, for the adjoint operator B*, 
the opposite characterization of all the eigenfunctions is in use: all the eigenf unctions 

(4.13) i^siy) of B* are not oscillatory and are of a "minimal" growth as y — )■ 00 . 

The last issue of ip%{y) being of a "minimal" growth (since there are many faster other 
asymptotics that are oscillatory) will be key in the nonlinear setting for the QLSE (I1.16p . 
In the linear case n = 0, all those notions admit a natural (but still not that easy) standard 
treatment, so we do not need to stress upon such issues in what follows. 

Thus, the definition (I4.12p of the operator B justifies the necessary part of the point 
spectrum with A/3 = — ^ and excludes its "positive part" (obviously nonexistent in (I4.12p ) 

(4.14) A+ = ^, |/3|>0, with ^;+ = ^/3. 

Nevertheless, for convenience, we will continue to refer to (I4.13P as a simple, efficient, 
and actually true way for a correct characterization of necessary eigenfunctions. In other 
words, the definition of B* actually includes special "radiation-like" conditions at infinity, 
which delete the non-desirable positive spectrum. 

4.4. Discrete spectrum and Hermitian polynomial eigenfunctions. Thus, defin- 
ing B* by (I4.12p . with the discrete spectrum only: 

(4.15) a(B*) = a(B) = {A/3 = -g, |/3| = 0, 1, 2, ...}, 

where all eigenvalues have finite multiplicity, and polynomial extended eigenfunctions to 
be determined explicitly shortly. Next, as customary, we fix other properties of the adjoint 
operator in a manner similar to Lemma 13.11 



Again, there are extended ones, ^ H^V^, and will restore a usual meaning in Appendix D. 
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Lemma 4.1. Under the above hypothesis and conditions: 

(i) (Extended) eigenfunctions ip^iy) are polynomials of order |/3| > 0; 

(ii) The subset of eigenfunctions $* = {ip^} is complete and closed in L^(M^); and 

(iii) $* is L'^,- evolutionary closed in the sense that the eigenfunction expansion l \4.11\) , 



which converges in the mean and uniformly on compact subsets, presents the rescaled 
{according to i \4-6^ ) solution of the LSE fli.ip for any data Uq G L^, (M^); and 

Proof, (i) Construction of polynomial eigenfunctions. Of course, the necessary discrete 
spectrum f l4.15p follows from (14. lip . We now intend to show how to obtain these results 
directly from the differential operator. 

Thus, ip*{y) is a polynomial. If its degree is k, then 

riy) = tPM^ 

j=0 

where Pj{y) is a homogeneous polynomial of degree k — 2m j with s = [^] , denoting the 
integer part. Since by the Euler identity 

N 

-hT.y^'-^ = -i,Po{y) = XPo{y\ 

we see that ^ = ~^ and Poiy) may be an arbitrary homogeneous polynomial of degree 
k. Other polynomials Pj{y) are then defined as follows: 



p,(y) = i(i(-Ar)^Po(?/), J 



s. 



This structure of 4'*{y) implies the completeness of the set of eigenfunctions in L"]^^^^ 



p^ 



In the second-order case m = 1, this construction leads to the generalized Hermite poly- 
nomials, which were already introduced in f l3.39p . Note that the polynomial structure of 
adjoint eigenfunctions follows from the expansion fl3.28p . where the coefficients of initial 
data Mp{wQ) in (IX^ are the dual products of Wq G Lp*(M^) and ^^ G ^^(M^). This 
implies that each ijj%{y) is a finite linear combination of elementary polynomials y'^ . 

We now fix Po{y) = y^ , so that, for (extended) eigenfunctions {'0/3} of B in (I3.34p . the 
corresponding adjoint eigenfunctions take the form 




(4.16) 

We also call (I4.16P the generalized Hermite polynomials. For m = 1, up to normalization 
constants, these coincide with those given by the classic generating formula (I3.39p . 

(ii) Completeness and closure. This is the well-known fact that polynomials {y^}, which 
are higher-order terms in any eigenfunction ?/)o, are complete in suitable weighted L^- 
spaces; see [511 P- 431]. Closure is associated with the eigenfunction expansion (14. lip , 
(iii) Evolution closure follows from (14. lip . D 
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5. Application I: Evolution completeness of $ in L^.(M^), sharp estimates 

IN ]R^"^\ AND SOME EXTENSIONS 

5.1. Linear PDEs. Our first result is about the following refined asymptotic scattering: 

Theorem 5.1. Consider the Cauchy problem f li.ip for uq G L^*(M^) and Uq ^ 0. Then 
there exists a finite I > and a function fi{y), such that, as t —> +cxd, 

(5.1) uix,t) = r^ h(^) + 0(r^)] 

uniformly on compact sets in y = -j^, where fi{y) is a nontrivial superposition of 
extended eigenfunctions {-0^, |/3| = /} of M from the corresponding finite- dimensional 
eigenspace. 

Of course, this is a corollary of our convergence analysis of the series fl3.20p and fl3.28p . 
where / is the minimal multiindex length \(3\ = I, for which Mg(i(;o) 7^ 0. As in Agmon's 
classic results for the parabolic case (see zero set applications of advanced Agmon-Ogawa 
estimates in [H] for parabolic PDEs for m = 1), a super-fast decay in (15.11) corresponding 
to / = 00 implies u{x, t) = 0, so that (other topologies are meant) 

(5.2) |w(x,t)| < t"-^ as t -)■ 00 for any i^ > 1 =^ u = 0. 

Further extensions of the above classification of t — )■ cxd behaviour are rather straight- 
forward for asymptotically small perturbations of the LSE such as 

(5.3) Mi = -i(-A)™M+ Yl a^ix,t)D^u in M^ x R+, 

(0<|7|<2m) 

with, say, bounded complex-valued coefficients {a^}, which decay sufficiently fast: 



(5.4) |a^(x,t)| = o(t2m 1) as t — )• -Fcx) for any < I7I < 2m — 1. 

We then deal with solutions u{-,t) G L^(R^), for which eigenfunction expansions make 
sense. Then after scaling (I3.13p . we arrive at the same equation (I3.10p with asymptoti- 
cally (and exponentially if necessary) small perturbations, which can be tackled by the 
eigenfunction expansion techniques, though some parts of the study can be indeed techni- 
cal. Some extra efforts are necessary to tackle convergence properties of such series that 
look rather technical, though can be involved in some places. 

5.2. NLSE: on a "centre subspace" behaviour. There are no doubts that, at least 
partially, the classification results for decaying as t — )■ 00 solutions can be extended to 
the NLSE (I1.15p . Indeed, the nonlinear term then has the form 

(5.5) i\u\^^'^u = ao{x,t)u, where ao(a:, t) = i |ii(x, t)|^"^ — )■ as t — ;■ +00 

on such small solutions. Therefore, one can expect that the perturbation techniques are 
also effective here. 

Agreeing with that and do not performing this routine, but sometimes technical anal- 
ysis, we would like to stress our attention to the principal fact showing that even this 
looking purely perturbation approach is not straightforward. Namely, we next formally 
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show that the NLSE (11.151) can admit some small solutions with a complicated asymp- 
totics corresponding to centre subspaces of the rescaled operators. 

Bearing in mind the typical linear behaviour (15. ip . where / = 0, 1, 2, ... , we perform in 
(11.151) the standard rescaling as in (13.131) . 

(5.6) u{x,t) = {l + t)-^v{y,T), y = jj^^^^j^, T = ln(l+t). 
The rescaled solution v{y,T) satisfies the perturbed equation 

(5.7) v^={B + ^l)v + ie-^'^\v\P-^v, where 7, = (£^^^±0 _ i. 

It follows that there exists a sequence of critical exponents {pi, I = 0,1,2, ...} such that 

(5.8) 7; = iff p = pi = l + ^^, l>0. 
In these critical cases, (15. 7p yields the autonomous equation 

(5.9) Vr = {B + ^I)v + i\vr^v. 

Then, we are looking for solution v{-,t) with the behaviour for r ^ 1 close to the 
centre subspace of the linearized operator B + ^ I. Such a centre subspace asymptotic 
dominance assumes that in the eigenfunction expansion of the solution 

(5-10) v{y, r) = J^^/s) C/3(r)^,3(y) 

the leading term for r ^ 1 corresponds to an eigenfunction Lpi of B, 

i.e., tpi belongs to the centre subspace of the linearized operator B+ ^^ / in the perturbed 
equation (15. 9p . Hence, we suppose that 

(5.11) v{t) = ai{T)[ipi + o{l)] as r -)■ +00 

and this asymptotic equality can be differentiated in y and r. Then the equation for 
the leading expansion coefficient ai{T) takes the form (recall the convention (14. 4p for the 
metric) 

(5.12) ai = i\ai\P'-^ai[ci + o{l)] for r > 1, where ci = {\^i\P'-^^i,^*), 

and Q 7^ as should be assumed. Here, for the first (and the last) time, we need the values 
of generalized extended linear functional to be introduced in Appendix D. We must admit 
that it is not easy to evaluate such coefficients q, even numerically in ID. Moreover, 
we still do not know any efficient method to get these values, so we will need further 
speculations. 

However, to confirm that such a formal analysis actually makes sense, we present a 
simple explicit example of such solutions on a centre subspace (in fact, on a manifold). 
Example: explicit centre subspace periodic solutions for m = 1, I = 0. Obviously, 
the simpler case is / = and m = 1, where the eigenspace is ID, so that (p^ = F and 
hence, by (I2.6p . 

2m 

V ^_ — /I ??IP0-1 T? J,*\ — I 

(47r)^/2 



(5.13) \F{y)\ = b, = -^, =^ c,= {\Fr-'F,ro) = b^' 
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so that Co is real. We have used the convention J F = 1, which, for A^ = 1, holds in 
the usual sense, since the improper integral converges. Then substituting into (15. 9p . with 
/ = 0, yields the following explicit solution: 
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(5.14) v{y,T) = ao{T)F{y) =^ do = i|aokoo =^ ^ |ao(r)p = 0, 

i.e., this explicit centre subspace behaviour is a periodic orbit. Note that this is not an 
L^(M^) or any L^(R^) solution, since F{y) = ip){y) is the first (extended) eigenfunction 
of the operator B, and also a standard eigenfunction of B in the extended space of eigen- 
function expansion closures; see Appendix D. Therefore, we do not know whether a centre 
subspace solution (I5.14p can have a stable orbit connection with more customary L^ and 
other solutions of Schrodinger equation theory. 

For general m > 2 and arbitrary / > 0, we need further arguments concerning the 
coefficients q in (I5.12p . For instance, using the continuity argument with respect to p, we 
again observe "almost" real values, since: 

(5.15) Q = 1 by orthonormality, ii pi = 1, so q ~ 1 if p; ~ 1 for all / ^ 1. 

However, such an asymptotic result assumes a technical proof, which falls out of the scope 
of the present analysis. Overall, for real coefficients q, the system (I5.12p is Hamiltonian, 

so, as above for / = 0, it describes a periodic behaviour close to the centre subspace. 

For complex valued q, such a centre subspace behaviour can be more complicated. 
For instance, for q = Iq, with q > (a rather hypothetical situation to be used as an 
illustration only), integrating (I5.12p yields the following decaying functions: 

(5.16) ai{T) = Cit'^{1 + o{1)) as r ^ +cx), 

where Q G C is a constant. In terms of the original (x, t, «)- variables, such a behaviour 
takes a form of a logarithmically perturbed linearized pattern 

(5.17) u{x,t) = Ciit\nt)-^[ipi{j^)+o{l)] ast^+oo. 

As we have mentioned, similar asymptotic patterns can be constructed for the "stable" , 
defocusing NLSE 

(5.18) Ut = -i{-/\Tu-i\u\P-^u in R^ x R (p > 1); 

see [m [HHl [92l [95] for key references and results concerning (15.180 for m = 1, as well as 
recent papers [5l [6l [331 ES [121 [121 1101] (and references/short surveys therein) for m = 2, 
i.e., for the biharmonic nonlinear Schrodinder equation 

(5.19) iut + /\^u = ±\u\P-^u in R^ x R. 

For (I5.18p . such a centre subspace approach looks even more promising than for the 
unstable PDF (I1.15P admitting also blow-up in these ranges. More flexibility is added 
when replacing the nonlinear term by a more general one, 

±i\u\^^~'^u ^^ d\u\^~^u^ where (i=a + i6GC, 067^0. 
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A rigorous justification of tlie centre manifold-like patterns of a periodic or fl5.17p -type 
is a difficult open problem, which we do not touch here. Notice that even existence of 
an invariant manifold (in which functional setting?-in an extended space of closures as in 
Appendix D?) is a very difficult problem for such Hermitian spectral theory dealing with 
the pair {B, B*}. This is regardless the good spectral properties of B listed in Lemma 
13.11 and also its sectorial setting in the topology of /^ in Proposition ID. 11 which however 
suggest a certain confidence that this behaviour can be verified by using the powerful 
machinery of classic invariant manifold theory, 



6. Applications II and III: Local structure of nodal sets and unique 

CONTINUATION 

6.1. Application II: blow-up formation of multiple zeros for linear PDEs (a 
Sturmian theory). Next, we arrive at the following classification of zeros of solutions 
oftheLSE flTTD : 

Theorem 6.1. Consider the Cauchy problem (1 1 . 1\] foruo G L^, (M^) anduo ^ 0. Assume 
that the corresponding solution u{x,t) creates a zero at a point (0,T), i.e., u{0,T) = 0. 
Then there exists a finite I > 1 and a generalized Hermite polynomial '^*p{y) such that 

(6.1) u{x,t) = {T-t)^[^^{j^^^)+oil)] as t^T- 

uniformly on compact sets in y = ,rp_^u/2m > where i~p*i{y) ^ is a superposition of poly- 
nomial (extended) eigenf unctions {f/^S, |/3| = /} o/B* from the corresponding eigenspace 
ker(B*-^). 

Since ipi{y) is a generalized Hermite polynomial, the multiple zero of Iieu{x,t) (or, 
equivalently, Imu^x, t)) occurs at the point (0, T~) by "blow-up focusing" of several zero- 
surfaces {x^{t)} of Re if* ( frp_^.i/2m ) ; which move according to the scaling blow-up law 

(6.2) x^(t) ~ (T-t)2^ ^0 as t -> T^ (ItI < 0- 

The result ( 16. ip follows from the series (14. lip , for which the rescaling (14. 6p is performed 
relative to the time moment T rather than 1. In a natural sense, the countable family 
of the types of asymptotics ( 16. ip describes the sharp "micro-turbulent" structure of the 
PDE (II. ip . since, by evolution completeness, on smaller space-time scales, the solution 
behaviour is trivial (a constant one mostly). In other words, ( 16. ip exhausts all possible 
micro configurations that can be created by the LSE (and also by many other related semi- 
linear and quasilinear PDEs admitting similar blow-up rescaling and Hermitian spectral 
properties; see below). See also |2l] for parabolic and other real- valued PDEs. 

6.2. Application III: Unique continuation. Various classic and other well-known new 
unique continuation results for linear and nonlinear Schrodinger-type PDEs can be found 
in [ini HH HJl So], where further references are available. These directions on uniqueness 
PDE theory have their origins in many principal works in the twentieth century including 
such a classic path as Holmgren (the starting point, 1901)-Carleman (1933)-Myshkis 
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(1948)-Plis (1954)- Calderon (1958)-Agmon-Nirenberg... . Here we present an example 
of a slightly different type of a "blow-up micro-scale uniqueness" study based on the 
spectral properties of B*, which is responsible for blow-up scaling of the PDE. 

Thus, all the types of nodal sets of zeros for the LSE (11 .ip are exhausted by the zero 
structures of Re or Im of all the generalized Hermite polynomials $* = {'0«(?/)} given by 
(14.161) including arbitrary linear combinations on all the eigenspaces. We fix this in the 
following rather unusual unique continuation theorem: 

Corollary 6.1. Let under the hypotheses of Theorem \6.1\ the nodal set of the real part 
{or, equally, of the imaginary part) of the solution 

(6.3) Af{u) = {xeR^, teR: Re u{x, t) = 0} 

has a nontrivial component that evolves as t ^ T^ in a manner that does not asymp- 
totically match the zero sets of any finite linear combinations of the real parts of the 
generalized Hermite polynomials from $*. Then 

(6.4) M(x,t)=0 m M^ X M. 

Of course, this is just a conventional version of the uniqueness result that is based on 
the eigenfunction expansion (14. lip , which can be expressed in a different and more reliable 
manner. For instance, e.g., depending on / and also much on m > 2, for the case, where 
the resulting polynomial Re Lp*i (y) does not change sign at all (so that (16.30 is locally 
empty), we then have to postulate just the asymptotic behaviour such as (16. ip . Then the 
alternative (and more correct and universal) sounding of the unique continuation result 
will be as follows: 

(6.5) if u{x, t) violates any of non-trivial asymptotics (16. ip near zero, then u = 0. 

Recall that, typically, for the real- valued evolution linear or nonlinear PDEs with in- 
terior regularity, unique continuation theorems stated in the pointwise sense deals with 
zeros of infinite order in the following manner: if, in a natural integral mean sense, 

(6.6) u has an infinite-order zero at (0,T), then u{x,t) = 0; 

see [TT] and [21] for further references and results for parabolic PDEs (such results are 
also known for the LSEs and are proved by iterating Carleman's classic estimates). Of 
course, (16. 6p becomes trivial for analytic solutions (though extensions to smooth non- 
analytic ones along the lines discussed below makes deep sense), so that we present a new 
pointwise uniqueness version (16. 3p . which looks not that trivial. 

It is natural to expect that the above classification of all the possible zeros remains 
for the perturbed LSEs such as (15. 3p . with arbitrary bounded coefficients {a^}. We then 
need to assume that u{x, t) is locally good enough close to the point (0, T~), and at least, 
u{-,t) e L^(R^) (see details in Appendix D), so we can use the corresponding eigen- 
function expansions endowed with a strong enough topology of convergence on compact 
subsets, which are sufficient to detect and identify the zero structure of solutions. At 
least, we need convergence a.e., which is guaranteed by the L^-metric. Nevertheless, the 
pointwise sense of such expansions will possibly demand some extra hypotheses that are 
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not discussed here. It is known that, even in the parabohc case, such a Sturmian theory 
on zero sets leads to a number of technical difficulties; see [21], where further references 
are given and applications to other classes of PDEs are discussed. Note also that such 
extensions can be applied to related partial differential inequalities (PDIs), e.g., for 

(6.7) \ut + i (-A)™?/| < C{\u\ + I Vm| + ... + l^^m-i^i)^ 

where C > is a constant. One can see that the right-hand side is always negligible after 
rescaling (14. 6p . so it does not affect the asymptotic zero classification (16. ip . 

6.3. The NLSE: similar local zero set behaviour. For the NLSE (ll.lSp . the local 
zero evolution remains unchanged since no centre subspace patterns are available. Namely, 
assuming again that a zero occurs at (0, T), we perform for (ll.lSp the scaling 

(6.8) u{x,t) = v{y,r), y = —H.^^, r = -ln(T-t) =^ v^ = B*v ±e-^ i\v\p-^v. 

In other words, the nonlinear term near the zero always creates an exponentially small 
as r — > +00 perturbation of the dynamical system for the expansion coefficients. Hence, 
it is very unlikely that this can somehow essentially affect the local zero structure near 
(0,T). Recall again that a rigorous analysis is rather involved even in simpler parabolic 
cases, [T71 [21]. Of course, close to the zeros, the NLSE f ll.lSp falls into the scope of the 

PDI dSZD- 
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Appendix A. "Schrodinger equations" are most popular in the twenty 

FIRST CENTURY ACCORDING TO MathSciNet 

We have claimed in Section [1] that Schrodinger PDEs are most popular currently in mathe- 
matical PDE theory. To "prove" this statement, we present below Table 1 on citations of the 
different types of equations in the MathSciNet. One can get these results by typing in the box 



Title I the corresponding type of equations, i.e., "hyperbolic (elliptic, parabolic, etc.) equation". 
Each resulting page in MathSciNet contains a variable number of papers (in average, it is about 
20, but can be more than that), so for convenience we count the number of pages. Naturally, at 
this moment, we put a blind eye on the obvious fact that there can be, and indeed many, papers 
on various PDEs, whose title do not contain those two words. Anyway, we believe that this is 
a reliable statistics and the table reflects some true information. 

Thus, the absolute winners for the total number of papers (via "All matches"), which have in 
the title "... equation" are "elliptic" and "parabolic" ones (recall that the MathSciNet operated 
that time with about 3.5 million papers, written by about 200 000 mathematicians). However, 
during two years 2007-08, more papers were published for "Schrodinger" ones: listing all of 
them takes about the record 25! pages. This adds extra bits to our motivation of developing a 
refined spectral theory for the Schrodinger PDE (jl.ip . 
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Table 1. Citations in MathSciNet, 8th September, 2008 

PDE type in the title All matches In 2007-08 

"hyperbolic equation" 3321 ~ 5 pages 

"elliptic equation" 7118 ~ 19 pages 

"parabolic equation" 6972 ~ 18 pages 

"Schrodinger equation" 5264 ~ 25! pages 

Appendix B. Application IV: Towards boundary point regularity theory 

We now present another application of a refined {B,B*}-spectral theory to the problem of 
regularity of boundary characteristic points for Schrodinger equations such as ()1.2p and (jl.ip . 
These regularity issues were always in the core of general potential theory, which in its turn 
represents one of the most classic directions of linear and nonlinear PDE theory initiated already 
by Dirichlet himself in 1820s. We refer to Maz'ya's monographs with collaborators |561 [57] for 
the history and key results on elliptic PDEs, as well as to recent surveys in |25[ [27] devoted 
more to parabolic PDEs, whose approaches and results will be essentially used later on. 

Of course, there are many very strong and classic boundary regularity results for Scrodinger 
equations (|1.2p and (|l.ip . which are explained in several key papers mentioned in the intro- 
duction. In particular, Schrodinger (or similar and often equivalent beam-type) equations in 
non- cylindrical domains have been studied in the 1960s by J.-L. Lions and E. Magenes, and by 
G.A. Pozzi, and later on by T. Gazenave and others. We refer to papers [3l[8l[l2j, where further 
references and results can be found. The existence-regularity results therein in principal cannot 
treat characteristic boundary points, where, as we show, the continuity of the solutions is rather 
tricky. We hope that our brief regularity exposition based on spectral theory of the pair {B, B*} 
will insert some new features to this classic area, which were not observed before. 

B.l. Regular boundary points and general asymptotic problem. Without loss of gen- 
erality, in order to explain the main ingredients of the boundary point regularity, we consider 
the simplest case m = N = 1, i.e., the ID second-order Schrodinger equation: 

(B.l) ut = iuxx in Qo, 

where Qq is a typical domain, for which (0,0) is a characteristic point, i.e., the straight line 
{t = 0} touches the lateral boundary of Qo- We define the backward parabola^ Qo as follows, 
assuming that it has a single finite right-hand lateral boundary (we will see why this is necessary): 

(B.2) Qo = {-oo <x < R{t), t£ (-1,0)}, 

where R{t) > for all i G [—1, 0) is a given sufficiently smooth function satisfying 

(B.3) i?(i)^0+ as t^Q-. 

Finally, we pose the Dirichlet boundary condition: 

(B.4) u{Rit),t)=0 for te(-l,0). 



In M.^ , it is a backward paraboloid with a quite similar study as for A^ = 1, though, of course, it becomes 
more involved; see those typical boundary regularity features in 28J to 3D Navier-Stokes equations. 
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and prescribe bounded smooth L^-initial data at i = — 1 (without loss of generaUty, we ahow 
uq{x) too decay exponentially fast as x — )• — oo or to even to be compactly supported): 

(B.5) u{x,-l) = uq{x) for x £ {-oo,R{-l)). 

Overall, ()B.ip - ()B.4p is a well-posed initial-boundary value problem, and we assume that it 
possesses a classic bounded solution up to the characteristic moment t = 0~ . 
Thus, the point (0, 0) is called regular in Wiener's classic sense if 

(B.6) u(0,0) = 

for any such data uq, i.e., there is continuity along the boundary, and irregular otherwise (if 
(IB.6|) fails for some uq). Our main goal is to show how to answer the following question: 



(B.7) for which lateral boundaries given by R{t), (0,0) is regular (irregular). 



In fact, this follows the canonical regularity statement by I.G. Petrovskii in 1934-35 |75ll76j . 
who almost completely solved the boundary regularity problem for the heat equation 

(B.8) ut = Uxx- 

This led to his famous "log log backward parabola", meaning, in particular, the following re- 
markable and delicate results: 



(B.9) 



(i) R{t) = 2 V^ v^ln|ln(-t)| =^ (0, 0) is regular, and 

(ii) R{t) = 2(1 + e)^/^ y^ln| ln(-t)|, e > =^ (0, 0) is irregular. 



In what follows, unlike many strong well-known results and approaches, we follow an asymptotic 
matching blow-up approach to the present regularity problem, which was developed in [25| and 
[27j for higher-order linear and nonlinear parabolic PDEs, respectively. 

Actually, of course, we are solving a more general problem on the asymptotic behaviour of 
solutions u{x, t) as a; — )• and t — )• 0~, so our goal is as follows: 



(B.IO) to describe sharp asymptotics of u{x,t) at the blow-up point (0,0 



and essentially detect their dependence on the function R{t) defining the right-hand lateral 
boundary. As usual, once (|B.10p has been solved, one can check the regularity property (|B.7p . 

B.2. A comment: t^vo boundary conditions must be non-Hamiltonian. Let us comment 
on the non-symmetric shape of the domain in ()B.2p . Using the symmetric shape (as can be done 
for the heat equation (jB.SP or for the bi-harmonic one (jl.Op ) 



(B.ll) Qo = {-R{t)<x<R{t), tG(-l,0)} 

is not possible. Indeed, since (|B.ip with the Dirichlet conditions on the lateral boundary of Qo 
is a Hamiltonian system with the L'^-conservation, as t — )• 0^ , one observes the concentration of 
this L^-energy onto shrinking-to-zero x-intervals, so that, obviously, 

(B.12) for domains (jB.lip and conservative boundary conditions, (0,0) is irregular, 

that eliminated the regularity issue at all. On the other hand, the symmetric domains such as 
(jB.lip are admitted if the boundary conditions on the lateral boundary violate the Hamiltonian 
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(symplectic) L^-conservation property. It is not that easy to find such conditions for the second- 
order equation (jB.ip . For instance, these could be the Robin ones (a third kind condition) at 
X = -R{t), 

Ux + cru = for some a £ C, 
and then the regularity would mean the continuity at (0, 0) along the right-hand boundary. 
For the fourth-order LSE, 

(B.13) ut = -iuxxxx, 

a domain (jB.lip also requires non-Hamiltonian lateral boundary conditions. By the identity 

(.^•J-4j ^ J Iti] — 11 ^UxxxU UxxxU^ _pj^r_i_\ + \UxxUx Uxx^x] _p^u\ 1; 

one can see that the homogeneous Dirichlet conditions 

u = ux = at x = ±R{t), t£{-l,0), 

are indeed Hamiltonian, the L^-norm of u[-,t) is preserved in such domains shrinking to a point, 
so (0, 0) is always irregular. The same is true for the Navier-type conditions 

u = Uxx = at X = ±R(t). 

Again, by ()B.14p . there exists the L^-conservation, so the vertex regularity problem makes no 
sense and (0, 0) is irregular for any nontrivial initial data. 
On the other hand, the following conditions 

(B.15) u = Uxxx = at X = ibi?(t), 

in general, violate the L^-conservation on shrinking domains a t — )■ 0~, so the characteristic 
boundary regularity problem makes sense. We can apply our "blow-up" scaling-matching ap- 
proach to study regularity of the vertex for conditions like in (JB.ISP or others of higher-order 
Robin-kind (but this will require a different boundary layer theory, see below). 

The boundary point regularity approach proposed here also covers the problem (IB.lSp . (IB.15P 
and other 2r?T,th-order ones with various (non-symplectic) boundary data. 

B.3. Introducing s\aw growing factor (/^(t). Thus, we return to the canonical (and indeed 
looking very simple) LSE-2 (jB.ip in the one-sided domain ()B.2p . with all conditions already 
specified. Then, similar to (|B.9p . we introduce a one-sided backward parabola at (0,0) given by 
the function 

(B.16) R{t) = y/^ lp{t), where r = - In(-t) -^ +oo as i ^ 0". 

In Petrovskii's criterion ()B.9p . 

c/?(r) ~ 2vlnr as r — t- +c«, 

so that ^p{t) is an unknown slow growing function satisfying 

(B.17) ifir) -^ -Foo, ip'{T) -^ 0, and ^ ^0 as r ^ -Foo. 

Moreover, as a sharper characterization of the above class of slow growing functions, we use the 
following criterion: 

(B.18) (^)'^oo as r^+oo {^'{t)^Q). 

This is a typical condition in blow-up analysis distinguishing classes of exponential (the limit 
is 0), power-like (a constant ^ 0), and slow-growing functions. See [831 PP- 390-400], where in 
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Lemma 1 on p. 400, extra properties of slow-growing functions (|B.18p are proved to be used 
later on. 

B.4. First kernel scaling. By (jB.16[) . we perform the similarity scaling 
(B.19) u{x,t) = v{y,T), where y = — ^. 

Then the rescaled function v{y,T) solves the rescaled equation 

'vr = 'B*v = ivyy - ^yvy in Qq = {-oo < y < (p{t), r > 0}, 

(B.20) <v = at y = c^(r), r > 0, 

^v{0,y) =vo{y) = uo{y) on (-oo,i?(-l)), 

where, by obvious (blow-up micro-scale) reasons, the rescaled differential expression B*, defined 
as in (jl.4p for ?tt, = A^ = 1, appears. In view of the assumed divergence ()B.17p . it follows that our 
final analysis will essentially depend on the spectral properties of the corresponding restricted 
linear operator B* on the whole line y G M, i.e., we arrive at the necessity of Hermitian spectral 
theory developed above. Note that, in particular, this differs our regularity analysis from several 
well-known ones such as Kondrat'ev's classic results of the 1966-67 [52l HS] (see a later survey 
|54]). where, as a rule, the rescaled boundary remains asymptotically fixed, which is key for 
using the spectral properties of the bundles of linear operators in locally compact domains. We 
will present further more detailed comments on that below. 

B.5. Regularity of a fixed backward parabolae is not obvious. First of all, we need to 
comment on the regularity of the vertex (0, 0) of the backward fundamental parabolae: 

(B.21) R{t) = l\/^, i.e., ip{T) = I = const. > 0. 

Then the problem (|B.20P is considered on the fixed unbounded interval 

(B.22) h = {-oo <y<l}, 

so that the final conclusion entirely depends on spectral properties of B* in Ii with Dirichlet 
boundary conditions. Since we need a sharp bound on the first eigenvalue, the clear conclusion 
on regilarity/irregularity becomes rather involved, where numerics are necessary to fix final 
details. In addition, as we pointed out, in more general setting for the fundamental backward 
paraboloids in IR^, the existence, uniqueness, and regularity of solutions in Sobolev spaces was 
proved in a number of papers such as [Ml EQI [Til [21], etc. Note that in [TQl p. 45], the zero 
boundary data were understood in the mean sense (i.e., in the L^-sense along a sequence of 
smooth internal contours "converging" to the boundary). 

Note that it is not obvious at all that the spectrum of B* on Ii (with a standard L^-setting as 
y — )• — oo) is real. However, one can expect that, by a continuity argument, the "first" eigenvalue 
-^0 = -^0(0 depending on / satisfies 

(B.23) Ao(/) ^0 as A ^ -Foo. 

This reminds a standard asymptotic fact from classic perturbation theory of linear operators 
(see Kato [46]) that the spectrum of B* in // approaches as I — )• -|-oo that in Lp(R), according to 
Lemma [4. 11 Then along with (IB.23p . one can also expect convergence of the first eigenfunction: 

(B.24) v*(y;0 ^ V'S(y) = 1 as i^+oo. 
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Notice that (|B.24p contains some features of a boundary layer that occurs as / — t- +cxd, which we 
will be use in the non-stationary limit (^(r) — )• +00 as r — )• +00. 

Overall, the limit ()B.23P reflects the possibility for (0, 0) to be regular or irregular for different 
values of / > depending of the sign of Re Ao(/) (or Xo{l) itself, provided that it is real). In other 
words, we conclude as follows: if the limit ()B.23P is oscillatory, then the backward parabolae 
()B.2ip can be regular or irregular. Note that this happens for the bi-harmonic equation (|1.9p . 
where the corresponding fixed parabolae with 

R{t) = {-t)k 

is regular for I = 4 but is irregular for / = 5. These conclusions for (jl.Op were fully justified 
numerically only, [251 § 6]. 

Thus, the regularity analysis of the backward parabolae ()B.2ip for the operator (jB.ip remains 
open, and its complete solution cannot be done without using enhanced numerical methods. 
Nevertheless, despite such a theoretical gap for constant Ts, we proceed to study the regularity 
for unbounded functions ()B.17p . which promises even greater mathematical challenge. 

B.6. The case (/?(t) — t- as r — )■ +00 is always regular. Indeed, in this case, as follows from 
(|B.20p . (0,0) is always regular, since (^(r) — )• and hence, by just the continuity of the solution 
v{y, r), we have v{0, r) — )■ as r — t- +00. It is worth mentioning that this is a completely rigorous 
result. Therefore, according to our blow-up asymptotic approach, the regularity problem of the 
vertex (0, 0) is trivial. Note that it is not that trivial via general PDE approaches (not including 
blow-up scalings); cf. much weaker assumptions in, e.g., |12j . 



B.7. Second scaling: Boundary Layer (BL) structure. Meantime, we return to the case of 
unbounded functions (/3(r)'s. Then, using standard boundary layer concepts of Prandtl-Blasius 
developed in 1904-08, we observe that, sufficiently close to the right-hand lateral boundary of 
Qq, it is natural to introduces the variables 

(B.25) ^ = 17^ and v{y,T) = w{z,t) =^ Wr = -^iW;,;, - ^ ZW;, + ^ zw^. 



ifiir) """" ""Vy ' ^ ~ "'V"'' ' '' ' ^'^ ~ 1^2 " ""^^ 2 "'""^ ' ip 

We next introduce the standard BL-variables (the same as for the heat equation) 

(B.26) C = ip^{T){l- z) = ip{y-ip), ip^{T)dT = ds, and w{z,t) = p{s)g{^,s), 

where p{s) is an unknown slow decaying (in the same natural sense, associated with (JB.ISP ) 
time-factor depending on the function fir). We will use later on the fact that 

(B.27) p = ao{T){l + o{l)) as t -^ +00, 

where ao(T) is the first Fourier coefficient of the solution v{y,T) relative to the adjoint basis 
i'ipl} of the operator B*. 

On substitution into the PDE in ()B.25p . we obtain the following perturbed equation (see 
details of a similar derivation in [251 § 7.2]): 

(B.28) g^ = Ag-jj{^ + ^)Cg^-fg^, where Ag = ig" + ^g'. 

As usual in boundary layer theory, we are looking for a generic pattern of the behaviour described 
by (JB.28P on compact subsets near the lateral boundary, 

(B.29) \^\=o{ip-'^{T)) =^ \z-l\ = o{ip~'^{T)) as t ^ +00. 
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On these space-time compact subsets, the second term on the right-hand side of (|B.28p becomes 
asymptotically small, while all the others are much smaller in view of the slow growth/decay 
assumptions such as (JB.ISP for (/?(r) and p{s). 

Then posing the asymptotic behaviour at infinity: g{^,T) is bounded as ^ — t- -|-oo, and 

(B.30) supg |(?(^,s)| = 1 (Hypothesis I for generic patterns). 

This is a typical "normalization by 1" condition from boundary layer theory. Note that, in view 
of highly oscillatory nature of any solutions of Schrodinger equations, we cannot normalize by 
1 any component of g{S,, s) = /i(^, s) + 'i'w{^, s). The condition (|B.30p will be used for matching 
with the solution asymptotics in the Inner Region. 

Then, as ^ — )• -|-oo, all the derivatives are assumed to be bounded, and we arrive at a standard 
stabilization issue of passing to the limit as s — t- -|-oo in ()B.28p . ()B.30p . Assuming that, by the 
definition in ()B.26p . the rescaled orbit {g{s), s > 0} is uniformly bounded, by classic parabolic 
theory [T8|, one can pass to the limit in (|B.28P along a subsequence {s/.} -^ -l-oo. Namely, by 
the above, we have that, uniformly on compact subsets defined in (|B.29p . as fc — t- oo, 

(B.31) g{sk + s) ^ h{s), where hg = Ah, h = at ^^ = 0, 

and (cf. (|B.30P ) h is bounded (i.e., being oscillatory) as .^ — )• -|-oo. The limit (at s = -|-cxd) 
equation obtained from (JB.28P : 

(B.32) hs = Ah = ihi:^ + ^h^ 

is a standard linear PDE in the unbounded domain M_|_, though it is governed by a non self-adjoint 
operator A. We then need the following its property: in an appropriate weighted L^-space if 
necessary and under the hypothesis ()B.30p . the stabilization holds, i.e., the w-limit set of {h{s)} 
consists of equilibria: as s — )• -|-oo, 

(h{^, s) -^ go{^), where Ago = for ^ > 0, 

l5'o = for ^ = 0, supg |5fo(^)| = 1 (bounded in M+). 

The stationary problem in (]B.33P can be easily solved to give the BL profile 

(B.34) ^o(^) = i(i_ei|), where |5o(0l = I cos(|)| < 1. 

We must admit that ()B.33p and (|B.34p actually define the class of solutions we are going 
to treat later on. Hopefully, this should be a generic class. So, we will not concentrate on 
the stabilization problem ()B.33p , which reduces to a standard spectral study of A in a weighted 
space. Actually, the convergence ()B.3ip and (IB.33P for the perturbed dynamical system (1B.28P is 
the main Hypothesis (H), which characterizes the class of generic patterns under consideration, 
and then (IB.30p is its partial consequence. Note that the uniform stability of the stationary 
point go in the limit autonomous system (IB.32P in a suitable metric will guarantee that the 
asymptotically small perturbations do not affect the omega-limit set; see p9l Ch. 1]. Such a 
definition of generic patterns looks rather non-constructive, which, however, is unavoidable for 
higher-order PDEs without positivity and order-preserving issues. 

B.8. Inner Region expansion: to^vards regularity. We next proceed as in \25\ § 7.3,7.7]. 
Namely, in Inner Region, we deal with the original rescaled problem ()B.20p . 
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In order to apply the standard eigenfunction expansion techniques by using the orthonormal 
set of polynomial eigenfunctions of B* given in (j4.16p . as customary in classic PDE and potential 
theory, we extend v{y, r) by for y > (/^(r) by setting: 

v{y,T) for < 2/ < c^(r), 



(B.35) viy,T) = v{y,T)Hi^iT)-y)-, ^ ^^^ y>^{T), 

where H is the Heaviside function. Since ?; = on the lateral boundary {y = 9?(r)}, in the sense 
of distributions, 

(B.36) Vr = VrH, Vy=VyH, Vyy=VyyH -Vy\^^^5{y -ip). 

Therefore, v satisfies the following equation: 

(B.37) Vr = ^*v + Vy\ ^ 6{y-ip) in R+ x M+. 

Since the extended solution orbit ()B.35P is assumed to be uniformly bounded in LL(M), we use 
the converging in the mean (and uniformly on compact subsets in y) the eigenfunction expansion 
via the generalized Hermite polynomials M.lGh : 

(B.38) v{y, t) = E(fc>o) MT)tPliy)- 

Obviously, this assumes the inclusion v{-,t) G L^, for all r » 1, which, by classic regularity 
theory for Schrodinger equations, is not a restrictive assumption at all. Then, substituting 
(IB.38|) into (IB.37P and using the orthonormality property (ID.7P in a "?;.p." sense of (|D.6p yields 
the following dynamical system for the expansion coefficients: 

(B.39) a'^ = Xkak + Vy\y^^^^^{6{y -ip{T)),tl;k) for all fc = 0, 1, 2, ... , 

where A^ = — 2 are real eigenvalues (J3.33P of B*. Recall that Afc < for all fc > 1. More 
importantly, the corresponding eigenfunctions tpkiu) are unbounded and not monotone for k > 1 
according to (14.161) . Therefore, regardless proper asymptotics given by (IB.39p . these inner 
patterns cannot be matched with the BL-behaviour such as (|B.30p . and demand other matching 
theory. Since these are not generic, the latter will be dropped, though can be taken into account 
for a full classification of (non-generic) asymptotics. 

Thus, one needs to concentrate on the "maximal" first Fourier generic pattern associated with 

(B.40) k = 0: Ao = and ^S(y) = 1 (V'o(y) = F(y)), 

where F{y) is the "Gaussian" profile (j2.6p . Actually, this corresponds to a naturally understood 
"centre subspace behaviour" for the equation ()B.39P : 

(B.41) v{y,T) = ao{T) ■ 1 + w {y,T), where w G Span{-i/;^, k>l}, 

and w-^{y,T) is then negligible relative to ao (''")• This is another characterization of our class of 
generic patterns. Hypothesis II. The equation for ao(T) then takes the form: 

(B.42) < = Vy\y=^^r)M^(^))- 

We now return to BL theory established the boundary behaviour ()B.26P for r S> 1, which for 
convenience we state again: in the rescaled sense, on the given compact subsets, 

(B.43) v{y,T) = p{s)go{ipHr){l - ^)) + ••. ■ 

Of course, since the limit BL-profile (]B.34p is uniformly bounded but essentially oscillatory, we 
are talking about matching of this BL-asymptotics with the constant one in ()B.4ip only in a 
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natural "average sense". E.g., after a standard integral averaging the oscillatory BL-profile, 
which eliminates the non-essential multiplier ^ in (IB.34[) . sincqj 

-j; /p cos^(|)d.^ — ;> ^ as L -> +00. 

Note that there are (finitely) oscillatory generalized Hermite polynomials V'fc(y) — Hk{iy) for 
large k (see ()4.16p . /c = /3 in ID), but these all are unbounded as y — )■ cxd, so cannot be matched 
with uniformly bounded BL-expansions. 

By the matching of both Regions, one concludes that, for such generic patterns, 

(B.44) ^^1 ^s T^+00 =^ p{s) = ao{T){l + o{l)). 

Then the convergence ()B.3ip . which by a standard regularity is assumed to be also true for the 
spatial derivatives, yields, in the natural rescaled sense, that, as r — )• +00, 

(B.45) Vy\y^^^^^-^ p{s)ip{T)-fi-^ao{T)ip{T)ji, where 71 = -5^(0) = §. 

We again recall that such an estimate is assumed to be true for a fixed above generic class of 
solutions satisfying a proper stabilization property in the BL. 

Thus, this leads to an asymptotic ODE for the first expansion coefficient for generic patterns: 

(B.46) a'o = 7iaoc^(T)e'^^ + ..., where 71 = ^^ = ^. 

This not-that-easy asymptotic ODE gives insight into main difficulties that one can face while 
posing and studying the problem on the boundary regularity of the vertex (0, 0) for Schrodinger- 
type operators. To this end, we first derive the real form of this system for 

(B.47) ao(r) = 6o(T)+ido(r), 

where these parts now satisfy the system: 

^0' = ^ [(^0 + do) cos (^) + (60 - do) sin (^)] + ... , 



(B.48) 



4' = ^[(-6o + do)cos(^) + (6o + do)sin(^)]+.... 

Note that the regularity of (0, 0) assumes that both limits are zero: 

(B.49) bo{T) and do{T) — )> as r — > -|-oo. 

The system (|B.48p shows a general range of problems that appear determining the conditions 
on the lateral boundary given by the function (p{t) to ensure (|B.49p . Of course, there is no any 
hope to guarantee (|B.49|) via a kind of Osgood (1898)-Dini-like integral condition of Petrovskii's 
type for the heat equation (|B.8p . The latter one has a very simple form (see a derivation in the 
same lines of spectral properties of B* in [251 § 7.7]): 



(B.50) 



(0,0) is regular iff J (p{s)e ^~ ds = -|-oo. 



which was already obtained by Petrovskii in 1934 [75 j [76] (earlier related results were due to 
Khinchin, 1924, in a probability representation; see [251 § 3.2] for further details). 

As a clue to some hard features of the system ()B.48p . consider as a "toy- model" a single 
equation of a similar form: 

(B.51) bo'^^Jocosi^). 



This is about averaging of I50P; averaging of |(7o| will give another (non-important) constant. 
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Integrating one obtains that 

(B.52) In \bo{T)\ ~ ^ / ^{s) cos {^) ds. 

In particular, the regularity of (0, 0) demands that 

(B.53) In |5o(t)| — )■ — oo =^ J (/?(s) cos ( "^ ^ ) ds — )■ — cxd as r — )• +c«. 

This immediately implies (this is true for the toy-model, but it seems such a condition exists for 
the whole system ()B.48P ) 

(B.54) J ip{s) cos C^ ^^' ) ds converges =^ (0)0) is irregular. 

It is easy to see that for the power-type functions 

(B.55) v{t) ~ t'^ =^ a > 1 means integral convergence and hence irregularity. 

Though fir) in ()B.55P is not a slow growing function, the above results are extended to that 
case. Hence, for any slow growing functions, the integral in ()B.54p always diverges. 

However, on the other hand, the divergence of the integral in (IB.54p . e.g., for a < 1 in 
()B.55p does not imply any regularity. Indeed, according to ()B.53p . the divergence must be 
to — oo, so that to create a regular boundary point (0,0) a special oscillatory cut-off of the 
boundary is necessary. This is such a "correction" of the shape of the boundary that eliminates 
the positive divergence part in the integrals in (IB.53p . Such a procedure, which establishes a 
coherent behaviour of the lateral boundary close to the characteristic point (0, 0) and oscillations 
of the rescaled kernel F{y), is inevitable for infinitely oscillatory fundamental solutions. This 
oscillatory cut-off already appears for the bi-harmonic equation (jl.9p : see [25l § 7.5], where 
further related details can be found. 

We hope that the above analysis correctly describes a full range of problems that appear in 
the regularity analysis and also correctly shows how {]B,B*}-spectral theory naturally enters this 
classic PDE area. 

Appendix C. Application V: Towards countable families of nonlinear 

EIGENFUNCTIONS OF THE QLSE 

C.l. The QLSE and its applications. This is a more striking and even controversial appli- 
cation of the spectral {B,B*} theory developed above. Namely, we now turn to a 2?7T-th-order 
quasilinear Schrodinger equation (the QLSE-2m) ()1.16p . Quasilinear Schrodinger-type models, 
some of which can be written as 

(C.l) iut + Au + f3\u\P~^u + e{A\u\'^)u = in M^ x M+, 

are not a novelty in several physical situations such as superfluid theory, dissipative quantum 
mechanics, and in turbulence theory. We refer to some papers from the 1970s and 1980s j38[ 
[79l 184] . to [34] for further references, and to |63] for a more mathematical knowledge, as well as 
to [64] and to Zakharov's et al papers [Ml [97], as a sufficient source of other reference and deep 
physical results. Note that the quasilinear model proposed in |64j in 1997 is more related to a 
kind of a "p-Laplacian operator" structure with fractional derivatives, such as 

(C.2) iui = A|L»,|^(||Z?,|^n|'|L»,|ln) + |Z?,rn, 
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with a standard Fourier-definition of operators l-D^I") fiaving the symbol |^|", so that, in par- 
ticular l-DxP = —D^ > 0. Then A = 1 in (1C.2P corresponds to the original defocusing model. 
Here, real parameters a and /? control dispersion and nonlinearity respectively. The standard 
NLSE then occurs for a = 2 and /3 = in ([02]) . 

For n = 0, the QLSE (J1.16p formally reads as the linear original one (jl.ip . Indeed, the QSLE 
()1.16p still remains rather rare and seems even an exotic equation in PDE theory, and therefore 
we now are going to propose a general approach to understanding of its new internal propertieso 

Overall, bearing in mind some clear discrepancies between the physical quasilinear models 
and the proposed one (I1.16p . we just say that this one was chosen as a typical example only 
to demonstrate our branching approach, so that others, more physically motivated models, 
would do the same, when the main mathematical ideas would have been properly explained and 
motivated. 

C.2. Two "adjoint" nonlinear eigenvalue problems. Thus, as in the linear case for n = 0, 
we are going to study global asymptotic behaviour (as t — )• -|-oo) and finite-time blow-up behaviour 
(as t^T- < -Foo) of solutions of the QLSE (fLlB]) . 

Overall, we are looking for similarity solutions of (I1.16P of two "forward" and Sturm's "back- 
ward" types: 

(i) global similarity patterns for t S> 1, and 

(ii) blow-up similarity ones with the finite-time behaviour as t — t- T~ < oo. 

Both classes of such particular solutions of the QSLE (I1.16P are written in the joint form as 
follows, by setting T = in (ii): 

(C.3) u±(x,t) = (±t)-"/(y), y = x/{±tf, where /3 = i^, for ± i > 0, 

where similarity profiles f{y) satisfy the following nonlinear eigenvalue problems, respectively. 



(C.4) (NEP)±: B±(a,/) = -i(-Ar(|/r/)±/3yV/±a/ = in 



aN 



Here, a G M is a parameter, which stands in both cases for admitted real (!) nonlinear eigenval- 
ues. Thus, the sign "-I-", i.e., t > 0, corresponds to global asymptotics as i — ?■ -|-oo, while "— " 
(t < 0) yields blow-up limits t — )• T = 0~ describing a "micro-scale" structures of the PDE. In 
fact, the blow-up patterns are assumed to describe the structures of "multiple zeros" of solutions 
of the QLSE. As we have mentioned, this idea goes back to Strum's analysis of solutions of the 
ID heat equation performed in 1836 [85j; see [23< Ch. 1] for the whole history and applications 
of these fundamental Sturm's ideas and two zero set Theorems. 

Being equipped with proper "boundary conditions at infinity" , namely, 

(C.5) for global case, B^(a,/) : f{y) is "maximally" oscillatory as y — )• oo, and 

(C.6) for blow-up case, B^(a,/) : f{y) has a "minimal growth" as y — )■ oo, 

equations ()C.4p produce true two nonlinear eigenvalue problems to study, which can be consid- 
ered as a pair of mutually "adjoint" ones. Note that ()C.6P actually also means that the admitted 



^While being in Bath in 2008, Peter Markowitz, answering a question of the first author, fast and witty 
called (I1.16p . m = 1, an NLSE of the porous medium type bearing in mind the classic parabolic PME 



ut = A(|u|"u) in M^ x E+, where n > 0. 
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nonlinear eigenfunctions are not of a type of a maximal oscillatory behaviour at infinity that 
connects us with the issue (J4.13P (which, however, is not sufficient, and a growth analysis at 
infinity should be involved in parallel, as shown below). 

Let us discuss in greater detail the meaning of those above conditions at infinity. Firstly, 
(jC.SP means that, due to the type of nonlinearity I/I", the oscillatory component such as (j2.1ip 
is admissible, with, of course, an extra generated algebraic factor of the WKBJ-type, which we 
do not specify hereby. This can be explained as follows: if f{y) has a standard WKBJ-type 
two-scale asymptotics 

(C.7) /(y) ~ |y| Vl^l" as y^oo, 

then, since |e"l^l"| = 1 for a G iM, substituting into (|C.4j) | yields the balance 

(C.8) Jn + (2m-l)(a-l) = 1, 

i.e., different from the purely linear one as in ()2.9p . since, for n > 0, the exponent 5 from 
the slower varying factor is involved (we do not calculate it here being a standard asymptotic 
procedure). 

Secondly, (1C.6P assumes actually also a "minimal" growth at infinity. Namely, quite similar 
to the linear problem for n = 0, the first two terms in (IC.4P generate a fast growing bundle: 
as y — 7- oo (as usual, we omit slower oscillatory components) 

(C.9) _i(_A)-(|/r/)-/3y.v/ + ...=0 =^ /(y)~|y|^. 

On the other hand, two linear terms in ()C.4p lead to a different slower growth as y — )■ oo: 

a 2 in I a I 

(CIO) ... -/3yV/-a/ = =^ f (y) ^ \y\~ ^ = \y\^^^ 

(recall that a(0) = A < 0). Since 

I'p ii^l 2m|a| 2m 

y^-^^) T+1^ < IT' 

this actually means that ()C.6P establishes a kind of a "minimal" growth of admissible non- 
linear eigenfunctions at infinity corresponding to (jC.lOp . For n = 0, this implies a polynomial 
growth, and all the admissible (extended) eigenfunctions of B* turned out to be generalized Her- 
mite polynomials (j4.16p . Note that, in self-similar approaches and ODE theory, such "minimal 
growth" conditions are known to define similarity solutions of the second kind, a term, which was 
introduced by Ya.B. Zel'dovich in 1956 [98], and many (but indeed easier) such ODE problems 
have been rigorously solved since that. For quasilinear problems such as ()C.4p . the condition 
()C.6P is incredibly more difficult. We thus cannot somehow rigorously justify that the problem 
()C.4p ■ ()C.6P is well posed and admits a countable family of solutions and nonlinear eigenvalues 
{a~{n)}. Actually, the homotopy deformation as n — )• 0^ is the only our original intention to 
avoid such a difficult "direct" mathematical study of this nonlinear blow-up eigenvalue problem. 

All related aspects and notions used above and remaining unclear will be properly discussed 
and specified. 

Of course, these conditions (1C.5J) and (IC.6J) remind us the "linear" ones associated with (I2.13P 
for B and (I4.13P for B* respectively, justified earlier for n = 0. Indeed, a better understanding 
of those conditions in the nonlinear case n > demands a much more difficult mathematics. 
However, one can observe that both (jC.SP and ()C.6P are just two asymptotic (not global ones) 
problems concerning admitted behaviour of solutions of ()C.4p as y — )• oo, so that, at this moment 
we are in a position to neglect these and to face more fundamental issues to be addressed below. 
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Note also that, at least, in ID or for radially syminetric solutions in R^, such asymptotic 
problems for not that hard nonlinear ODEs are easily solvable. 

Thus, for n = 0, equations ()C.4p . equipped with proper weighted L^ spaces, take the very 
familiar form: the corresponding differential expressions are 

(C.12) B, B* : -i(-A)™/±/3yV/±a/ = in M^. 

Then we observe the obvious relation between a's and A's from the spectrum a(M>) = a{M*): 

(C.13) for B : q = -A + ^ and for B* : a = X. 

Thus, our next goal is to show, by using any means, that, at least for small n > 0, the 
nonlinear eigenvalue problems 

(C.14) (NEP)-i- admit countable sets of solutions $^(n) = {a^, /^}|^|>o, 

where, as usual and as it used to be in the linear case, 7 is a multiindex in R to numerate the 
pairs. 

The last question to address is whether these sets 

(C.15) $ (n) of nonlinear eigenfunctions are evolutionary complete, 

i.e., describe all possible asymptotics as t — )• +00 and t — ;■ 0" (on the corresponding compact 
subsets in the variable y in ()C.3P ) in the CP for the QLSE ()1.16p with bounded integrable (and 
possibly compactly supported, - any assumption is allowed) initial data. Our main approach is 
the idea of a "homotopic deformation" of ()1.16p as n — )■ 0+ and reducing it to our linear equation 
(II. ip . for which both problems (IC.14P and (1C.15P are solved positively by a non-standard and 
not self-adjoint spectral theory of the linear operator pair {B,B*}. 

C.3. Example: first explicit nonlinear eigenfunctions. For m = 1, the problem ()C.4p , 
has the following first pair associated with the explicit kernel ()2.6p : 

(C.16) "^W = 2Tk and /+(y) = (^)^e^. 

Not that surprisingly, regardless the degeneracy of the QLSE (jl.lGp at the zero level {u = 0}, 
(|C.16p shows that the solution exhibits no finite interfaces. This is in a striking difference with, 
say, as a typical example, the TFE-4 

(C.17) ut = -V- (l/rVA/) in R^ x R+, where n > 0, 

which is known to admit compactly supported solutions in both the FBP and the Cauchy 
problem setting; see [301 [3T] and, respectively, [l9l[20], as a source of main results and further 
references. 

Concerning the "blow-up problem", for any m > 1, the adjoint nonlinear eigenvalue one 
()C.4p has the obvious first pair 

(C.18) aoW = and fo{y) = l, 

i.e., the same as for n = 0, where V'o(y) = Ij as the first Hermite polynomial; see (|4.16p . 

Anyway, in our further analysis, we cannot rely on any explicit representation of any nonlinear 
eigenfunctions. So, we now very briefly begin to explain our approach. 

42 



C.4. Branching for the "for^vard" problem for i ^ 1. We perform our bifurcation- 
branching analysis following the lines of classic theory [71 [HI [58l [90], etc., in the case of finite 
regularity. However, we must admit that this classic one does not cover rigorously the type of 
n-branching as n — t- 0, which requires extra difficult study, which we cannot address here. 

Thus, we are looking for a countable set of nonlinear pairs {a;|(n), /+(y; n)}|^|>o, which are 
assumed to describe all (hopefully, for n > small) asymptotic patterns for the QLSE (|1.16p 
as t — > +00, up to possible centre manifold patterns; see Section [5^21 However, as we have seen, 
such special patterns do occur, if there is a certain interaction (a transitional behaviour) between 
the linear and the nonlinear terms in equations such as (j5.7p . This looks rather unreal for the 
current problem under consideration. However, the evolutionary completeness of the nonlinear 
patterns $"'"(n) as t — >• +00 remains an extremely difficult open problem, which probably will 
be extremely hard to solve completely rigorously. 

Thus, assuming that n > is small, we perform asymptotic expansions in the operators and 
coefficients in ()C.4p , . Evidently, the crucial one is in the nonlinearity, which requires: 

(C.19) |/|"/^/e-i-l^l=/(l + nln|/|+o(n)) as n^0+. 

It is clear that the neighbourhoods of the nodal set of f{y) are key for ()C.19P to be valid in any 
weak sense (precisely this is needed for the equivalent analysis of the inverse integral compact 
operators involved, where a proper justification must take place). If f{y) has a nice nodal 
set consisting of a. a. isolated and "transversal" a.e. zero surfaces {or just some points only), 
with no "thin" concentration subsets, the expansion (1C.19P can be valid even in the standard 
pointwise sense, or at least in the weak sense. However, we do not know and currently cannot 
prove such deep properties of the nonlinear eigenfunctions involved. Note that, for n = 0, the 
generating formula (I3.34p , with a proper knowledge of such nice zero set properties of the rescaled 
fundamental kernel F{y) (this is doable), guarantees such necessary properties of eigenfunctions 

i'-yiy)- 

The above discussion establishes the main hypothesis to make our branching analysis to be 
(almost) rigorous, being applied, of course, to the equivalent integral equation, where establishing 
some further compact and other necessary properties of the nonlinear integral operators would 
take some time, indeed. This can be also done, since the problems ()C.4p _[_ can be reduced to a 
semilinear form. 

The rest of the expansions in (|C.4p | are straightforward (here, we already fix by 7 an n-branch 
we are going to trace out): 

a^(n) = a^(0) + a^n + ... = — A^ + tt^ + a^n + ... , 
(C.20) , _ ,„,„ , . , N 



l—O'y {n)n 1 



T5— ■ - a^n + ... , where a 



2m 



where we have already omitted all o(n)-terms (this again assumes extra regularity hypothesis 
already discussed above). In the first line in ()C.20p . the parameter d^ is an extra unknown. 
Substituting all the expansions into ()C.4p | yields the following perturbed problem: 

(C.21) {M-\^I)f + nh + ... = 0, where h= [i (-l)'"+iA'"(/ln |/|) - a^y • V/ + a^/]. 

Thus, as n — 7- 0+, we must look for a solution close to the eigenspace 

(C.22) ker(B-A^I) = Span{V'7, I7I = 0- 
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Therefore, under prescribed hypothesis, solutions take the form 
(C.23) f = T,\a\=iCaip<7 + n4>^ + ..., 

where the expansions coefficients {co-}|o-|=/ and the orthogonal part (j)^ are unknowns. 
Finally, substituting ([a23|) into (jU%]) yields the 0(n)-problem 

(C.24) (B - X^I)(j)^ + h = 0. 

Thus, the necessary (and, in properly regular cases, the sufficient) orthogonality condition of 
the solvability of ([gM]) . 

(C.25) /i±* ker(B-A^/), 

yields the following Lyapunov-Schmidt scalar branching equation: for all \S\ = I, 

(C.26) (i(-ir+iA™(Ec.^.ln|Ec.^a|)-a7yVEca^. + «^Ec.^^, ^^^=. = 0, 

where an integration by parts in the second term can be performed to simplify the expressions; 
see the next Appendix D for further details concerning the meaning of such extended generalized 
linear functionals. We must admit again that a proper well posed way for calculating values of 
such extended functionals is not available, so we present such a branching analysis just as an 
example reminding analogies with the classic approaches. 

Overall, ()C.26P is the required algebraic system for the unknowns a^ and {co-} (a convenient 
normalization condition on the latter expansion coefficients may be added). As usual, once the 
branching equation ()C.26P has been properly solved, this allows one to get the corresponding 
unique solution of the differential equation ()C.24p . etc. 

Indeed, ()C.26P is a very difficult algebraic system, which is not of any variational form, so 
one cannot use powerful category-genus theory [71 [58] to predict a number of solutions, i.e., 
a number of such n-branches originated from the given eigenspace in (IC.22p . Note that, for 
m = 1, when there exists some "symmetry" of differential forms B, B* reflected in (13. 8p and 
(I3.27p . the system (IC.26P reveals some "variational- like" features, since the eigenfunctions ip-y 
and polynomials i/j* can be identified in a weighted space. However, this is supposed to happen 
in a space with indefinite metric, so we do not check how this can be helpful. 

It is worth mentioning that the problem on a sharp estimate of a number of solution branches 
emanating from an eigenspace, remains essentially open even for classes of well-understood 
variational operators. On one hand, the first conclusion is classic: the number of branches is not 
less than the dimension of the eigenspace: indeed, since the corresponding algebraic system (like 
(IC.26P ) remains also variational, the category of the functional set is not less than the linear 
eigenspace dimension, whence the result. But, obviously, a sharper estimate of the solutions 
number becomes essentially nonlinearity-dependent, so that this is not (and, possibly, cannot be 
in the maximal generality) completely understood. For a number of branches that can emanate 
from the trivial solution, there have been obtained some specific examples only. There exist some 
results for potential operators (see [H] and ^80j as a guide), and a very few for non-gradient 
and non-self-adjoint operators [61]. For the TFE-4 (jC.lTp . such a branching analysis as n — )• 0"^ 
reveals a lot of technical difficulties, though some problems for simple and semisimple eigenvalues 
are shown to admit a rather definite path towards, 0. 
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C.5. Blow-up scaling problem: t — t- 0^. This is quite similar. Since the nonlinear eigen- 
value equations (|C.4[) _[_ differ by the sign in the linear terms only (but the boundary-radiation 
conditions at infinity are entirely different), instead of ()C.2ip . we arrive at 

(C.27) (B* -A^/)/ + n/i + ... = 0, where /i = [i (-l)"^+iA'^(/ln |/|) + a^y • V/ - a^/] . 

Therefore, looking for solutions ()C.23p . now over the kernel of the adjoint operator (B* — A^I), 

(C.28) / = Eh=/C-C + "</>; + -, 

we arrive at the "adjoint" algebraic system: for all \5\ = I, 

(C.29) (i (-l)-+iA'"(E c.V* In I E Carj) + a^y^E c^C - "7 E c^C i^s)* = 0, 

which are not easier that the first one ()C.26p . and does not have any variational structure, so 
the same principal difficulties on the solvability (this is easier to do in the ID or the radial case) 
occurs, and especially on the number of solutions. 
Recall that the patterns 

(C.30) ^(x,t)~e"-W-EH=/C<xC(y) + -, where y = ^-^, r = -ln(-t), 

are assumed to describe, for all finite multiindices 7 and all admitted solutions of the branching 
equation ()C.29P , the whole variety of "micro-scale patterns" , which are available for the QLSE 
(|1.16p at a given arbitrary point (0, 0~), unless some centre-subspace-type patterns might appear 
due to nonlinearities involved. In any case, we believe that ()C.30P describe (at least, a.a.) of 
generic formations of "multiple zeros" of solutions. In other words, for small n > 0, multiple 
zeros of the Reu{x,t) and/or Imu(x,t) as (x,f) — )• (0, 0~) are created by a self- focusing of 
zero surfaces of the corresponding profiles f{y) given in ()C.28p . which comprises a proper linear 
combination of the generalized Hermite polynomials (j4.16p . Indeed, this study requires further 
extensions, however, a complete and fully rigorous answer seems cannot be achieved. 

With such an "optimistic point", we end up the list of possible applications of our refined 
scattering linear spectral {]B,B*}-theory for 2mth-order rescaled Schrodinger operators. 

Appendix D. Eigenfunction expansions and little Hilbert spaces 

Given below further developing of spectral theory is not necessary for our main applications 
concerning classification of all the global and blow-up asymptotics, which follow from the corre- 
sponding spectral decompositions of semigroup representation of solutions. However, we think 
that these accompanying results are interesting and actually allow to extend the results to wider 
classes of solutions. 

Recall that the main difficulty with a proper definition of the operator B (13.30p by its spectral 
decomposition (13.29P was associated with the fact that the necessary for us eigenfunctions {ipp} 
were extended, i.e., did not belong to the present domain. In other words, these eigenfunctions, 
which inevitably appeared in eigenfunctions expansions (|3.28p of rather "good" solutions of 
Schrodinger equations, were much "worse" that the solutions themselves. To get rid of such 
a controversy and to restore the true meaning of the operator pair {B, B*} (instead of its 
restriction {B,B*}), we first introduce new extended spaces of closures. 
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D.l. Subspace where <l> is closed. Given the complete subset $ = {tpp} for the non self- 
adjoint operator B, we define the linear subspace L^ of eigenfunction expansions, 

(D.l) f G Z^ iff V = X^C/sV'^ with convergence in Lp(R^), 

as the closure of the subset of finite sums 

(D-2) {E|/3|<xC/3^/3,^>0} 

in the L^-norm. 

For clarifying such a space, we now derive better estimates to see which {cp} satisfy (jD.ip . 
Namely, we will use the equality ()3.17p for the key exponent in the exponential representation 
(|2.1ip of the asymptotics (which is responsible for the sharp estimate (|3.16p of eigenfunctions 
ipjj). Then similar to (|3.25p . but sharper, we then obtain for I = \/3\ ^ 1, 



/ , s _/ 2l(a-l) ^cyr . Al{a-L) 2I l(2-a)r 2m-l -, ' 

(D.3) !p\M'-{^ i^^[^^] - (2H'^ = r^^[ ,^(j^),,^,,,,J ^-^ 

Next, since 

(D.4) \\vf^2 = E C/3C^ / /5V'/3V'7' 

" (/3,7) M^ 

(jD.SP implies that: 

(D.5) v e L'j, if ci3 does not glow faster than Z^^^""^) for / = |/3| > 1, 

where e > is any arbitrarily small constant. Of course, ()D.3P defines more optimal and weaker 
inclusion conditions, but ()D.5P clearly explains how this works. 

D.2. Bi-orthonormality of the bases. Obviously, this is a principal issue for all the applica- 
tions, where eigenfunction expansion techniques take part. As we have seen, the eigenfunctions 
expansions such as (I3.32P introduce standard linear functionals {wq^i^V), which are well defined 
for all functions wq € L^, (M ), so that, as usual, V'« is an element of the adjoint space Lp(]R ), 
with yO = 4r, as customary. 

As the next step, according to our construction above, we have to define some generalized 
extended linear functionals from the adjoint space L^* . On one hand, this would correspond 
to a standard procedure of extension of such continuous uniformly convex functionals by the 
Hahn-Banach classic theorem in linear normed spaces, [51] • As we have seen earlier, those linear 



functionals are well defined according to ()3.22p in L?, with the standard (not in any v. p. or a 
canonical regularized, etc., sense) definition of the integrals. 

On the other hand, such extended linear functionals cannot be understood in a standard 
sense, so that we refer to them as to generalized ones. It seems, a full proper definition of such 
extended linear functionals in a usual functional framework will require a deeper analysis of the 
actual functional spaces and metric/topologies involved, which will essentially decline us from 
main PDE applications [Especially, since in some of the applications, we do not and even cannot 
pretend to be mathematically rigorous.] 

Therefore, in other words, for any v G Z^, we define extended linear functionals for any /3 as: 

(D.6) (vjipl)^ = (v, V'rt) denotes the expansion coefficient Cjs of v in (ID.ip . 
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In view of the performed construction of the space L^ via closure of finite sums (jD.2p , it is not 
difficult to see that such generalized continuous linear functionals are defined uniquely (in view 
of the density of finite sums (|D.2p ). 

Overall, in the sense of (ID.6p . the standard bi-orthonormality of the bases {'ip/3} and {ipt} 
becomes trivial: 

(D.7) (ipp^ip*)^ = {ipi3,'ipj) = (5/37 fo^ ^"^y f^ a^d 7, 

where (•, •) is the usual duality product in L^(]R ) and Sp-y is the Kronecker delta. 

Similarly, using the subset $* = {ip%} of the generalized Hermite polynomials (I4.16p . we are 

obliged to define the corresponding subspace L^^ of eigenfunction expansions, and eventually 
treat similarly the adjoint extended linear functionals {w,ipp)^ for any w S -^p,*- 

In writing ()D.7p . we use the standard L^-metric, which is convenient to see the normalization 
since ■i/'/j is essentially the D" derivative and ip*o is a polynomial, so that, including normalization 
factors yields 

(D.8) (#,^^) = ^^^^ / D^F{y) {yf^ + ...) dy = 1 for any multiindex /3 

in the sense of formal integration by parts. Therefore, the integral itself is not of a standard 
meaning, but can be treated in an involved "?;.p.-like sense" , which is difficult to clarify, and we 
do not feel any actual necessity to do this. A similar formalism exists for the whole eigenfunction 
set occurred in ()D.7p . Thus, we will use ()D.7p in the eigenfunction expansions to follow, bearing 
in mind its actual meaning specified above in (ID.6p . 



D.3. Little Hilbert and Sobolev spaces. It is convenient to introduce a little Hilbert space 
/p of functions v = Y^ api^p G ip 



ti of functions f = X] o-pi^p ^ -'^p with coefficients satisfying 



(D.9) i:\ap? < 



CXO, 



where the scalar product and the induced norm are given by 

(D.IO) {v,w)Q = Y,apCp, w = Y,Cfi'iljp ell and \\v\\1 = {v,v)q. 

Therefore, $ is now treated as a Riesz basis in Lp(M ), [9l[32]. We next define a little Sobolev 
space h?!^ of functions v £ l^ such that 

X]|A;3C;3|2 <00. 

The scalar product and the induced norm in /i^™ are 
(D.ll) {v,w)i = {v,w)o + {Bv,Bw)o, Ml = {v,v)i = Z {^ + \^p\^)\cp\^ 

where our bounded operator B : h^"^ — )• l^ has the meaning B : {cp} — ;■ {Xpcp}. This norm is 
equivalent to the graph norm induced by the positive operator (— B + al) with a > 0. Then 
/ip™ is the domain of B in l^, and, by Sobolev's embedding theorem, 

(D.12) /i2™ c Pp compactly, 

which follows from the criterion of compactness in P, [51]. 
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D.4. Basic properties in /^. Firstly, it follows that B is self-adjoint (symmetric) in l'^, 
(D.13) {Bv,w)o = {v,Bw)o f or sdl v , w £ hi"" . 

Secondly, we state some other straightforward consequences. 

Proposition D.l. (i) The Hilbert space /^ is a dense subspace of L^ in L?SM.^); 

(ii) <^ = {V^/^} is complete and closed in I'j, in the topology o/Lp(R^); 
(iii) the resolvent (B — A/)~"^ for A o'(B) is compact in I'i; and 
(iv) B is sectorial in P . 

Proof, (i) /p is separable and complete since the same is true for the isomorphic Hilbert space 
t^ of sequences. Let us show that Z^ C L^. For any v E /p, 

/pbl'dy = /|E«/3^/P#l'dy<EI«/3pE/e-l^l1;^i^My)rdy, 
and by the same estimates as in (j3.23p . (j3.24p . we conclude that, for / = |/3| » 1, 
(D.14) Je-ls'l"|^D'^F(y)|^dy</'(-^+^), 



2-a 



where v = -^ > and e > can be an arbitrarily small constant. Therefore, li C Lz,. 
Concerning the density of Z^, we note that given a t; = ^ a/3'0/9 £ Li% the sequence of truncations 
{YllBlKK ^I3'^l3^'^ = 1,2,...} C /p converges to v in the topology of Lp(R^) as iC — )• oo by 
completeness and closure of {ip/s}- 

(ii) Since ^ is orthonormal in /p, it follows that the only element orthogonal to {ip(3} is 0, 
and hence completeness of {ipp} in I'i follows from the Riesz-Fischer theorem. It is closed as an 
orthonormal subset in a separable Hilbert space [51j. 

(iii) For any v = Y^ apil)p E /p from the unit ball Ti in ^p with X^ |a/3p < 1, 

(B — \I)~^v = Y^ bj^ipis, where 
(D.15) &/3 = ^=-^^7i^ = -^[l + 0(^)] for|/5|»l. 

Therefore, for any e > 0, there exists K = K(e) > such that for any v £Ti, 

By the compactness criterion in P |51j . (B — A/)~^ maps Ti onto a compact subset in /p. 

(iv) Recall that (B — A/)~^ is a meromorphic function having a pole ~ -j^ as A — t- since 
Ao = has multiplicity one [32]. We then need an extra estimate on the resolvent, which is 
easy to get in I'i (one can check that it is not that easy in the big space L?). In the sector 
r^e = {A E C : A / 0, |arg A| < § + 6*} with a 6* E (0, f ), for any v = Yl aisi^p E if,, we apply 
(jD.lSp by using the fact that i^ ''_^| < -m- sin 9 in Qg to get 

||(B - Xiy^vWo = (E la^Pp^)^ < ii^A'l^llo- 
Since B is closed and densely defined, it is a sectorial operator in /p, see [22]. D 



Similarly to (jD.ip . for the adjoint operator B* we define the subspace L"^ C. -Lp(M^), where 
the eigenfunction subset <1>* is closed, ti^^, ^p™; etc., and next continue develop similar theory 
of self-adjointness and other properties 
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